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PREFACE 


Last year I gave a series of Extension Lec- 
tures on non-linear partial differential equations 
at the Lucknow University. These have now 
been written down in the form of a monograph 
to be published by the University. This small 
volume aims at giving a preliminary account of 
the methods employed and the results achieved 
for non-linear equations of the standard types. 
The first two chapters contain generalisations 
of my previous results, and are being published 
now for the first time. These lectures were 
originally advertised under the heading “The 
Theory of Non-Linear partial Differential Equa- 
tions” but the present title has been considered 
more appropriate for the monograph. 

I take this opportunity of offering my grateful 
thanks to Professor Birbal Sahni, F.R.S., for 
kindly arranging for the lectures in Lucknow, 
and to Dr. A. N. Singh, D.Sc., for his kind help 
during the printing of the book. 

Jme 21, 1939 Raziuddin Siddiqi 
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INTRODUCTION 


The theory of the boundary value problems 
in linear partial differential equations has been 
extensively studied, and the results incorporated 
in a large number of memoirs and treatises which 
have been reported in the Enzyklopadie articles* 
by Sommerfeld and Lichtenstein. 

The non-linear partial differential equations 
were, however, neglected for a long time, owing 
probably to the fact that not many physical pheno- 
mena were known which gave rise to boundary 
problems in such equations. But the situation 
is changing rapidly, and modern explanations of 
various physical phenomena are yielding more 
and more non-linear equations. This fact has 
induced several investigators to take up the subject 

* A. Sommerfeld ; “ Randwertaufgaben in der Theorie 

der partiellen Differentialgleichungem Enzykl. d. Math. 
Wiss. II A. yc, pp. 504-570. 

L. Lichtenstein: “Neuere Entwicklung der Theorie 
partieller Differentialgleichungen zweiter Ordnung vom 
elliptischen Typus”: Enzykl. d. Math. Wiss. II c. 12 
pp. 1277-1334. 
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systematically, and thus two diflFerent methods 
of attack have been evolved. 

The first method is that of Green’s function 
well known in the potential theory. This method 
was developed by M. Emile Picard in 1890, and 
applied by him to linear as well as non-linear 
elliptic and hyperbolic equations. 

In 1907 Holmgren and Levi simultaneously 
developed a similar method for the linear parabolic 
equations, employing, of course, the newly dis- 
covered theory of integral equations. Inspired 
by the works of Picard and Levi, M. Maurice 
Gevrey made a fairly complete investigation of 
the parabolic equations in 1915-1914. 

The method of Green’s function has been 
well developed and applied to non-linear equations 
in a large number of memoirs and treatises, a 
selection of which is given in the bibliography. 

Now in physical theories as well as in mathe- 
matical analysis, the solution of a boundary value 
problem in the form of a Fourier series or of a 
series of eigenfunctions is of a special value. In 
fact, since the classical work of Fourier himself 
the linear problems have been often solved in 
this way. Lichtenstein was probably the first 
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writer to give the solution of the non-linear hyper- 
bolic equation 

dx^ d/2 - " 

in the form of a fourier series* 


u{x, /) — 2 t>n {t) sin nx. 

n 

At his suggestion the present writer under- 
took a detailed study of the non-linear parabolic 
and hyperbolic equations from this point of view. 
A few results have been published in papers 
mentioned in the Bibliography. 

A fairly generalised theory is being presented 
in these lectures. The first chapter gives the so 
called “Fourier Method” for the non-linear para- 
bolic equation : 


d 

dx 




The solution of the first boundary problem is given 
in the form of a series of Sturm-Liouville eigen- 
functions, whose asymptotic expansion plays a 
fundamental part in the treatment. The Fourier 
coefficients are determined with the help of an 


* Unless otherwise stated, the summation is always to 
be taken from i to oo throughout this book. 
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infinite system of non-linear integral equations 
which is solved by the method of successive 
approximations . 

The method is then shown to hold for equations 
of higher order and for elliptic-parabolic equations. 
It is further shown that the problems for non- 
rectangular domains and for mixed boundary values 
can be solved by reducing them to the above 
problem with the help of suitable transformations 
of variables. Systems of differential equations can 
also be treated in this way. [Cf. Siddiqi, (7)]. 

The second chapter gives the Fourier method 
for non-linear hyperbolic equations. The method 
is similar to that for the parabolic equations, but 
there is a fundamental difference. The parabolic 
equations have a unique solution in Fourier series 
for any domain arbitrarily large, whereas the 
hyperbolic equations have a solution for a sufficiently 
small domain. This is compensated by the fact 
that the boundary values for parabolic equations 
are more restricted than those for the hyperbolic 
equations. The generalisation for non- vanishing 
boundary values is given in the last section, though 
other generalisations for non-rectangular domains, 
for mixed boundary problems, for equations of 
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higher order etc., have been left out, their treatment 
being completely analogous to that for the para- 
bolic equations. 

Many problems have still to be considered 
both for the parabolic and hyperbolic equations.* 
For instance, if the right hand member of the 
given differential equation contains the derivative 

then the system of integral equations which 

determines the Fourier coefficients is replaced by 
a system of integro-differential equations, for which 
a method of solution has yet to be developed. 
Similarly, the second boundary problem in which 

the derivative vanishes both at x = o and at 
dx 

X = K still remains to be considered. 

The account for non-linear elliptic equations 
in Chapter III is taken from the classical memoir 
of M. Picard published in the Journal de Mathe- 
matique in 1890. Similarly the account for non- 
linear parabolic equations in Chapter IV is taken 
from M. Gevrey’s paper published in the same 

^During the time the book was in the press, these 
problems have been dealt with by S. M. Sundaram, [i], [2], 
[3], working at Madras and Hyderabad. [See Chap. 
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Journal in 1913. It will be seen from these 
Chapters that the method of Green’s function 
establishes the existence of the solution for suffi- 
ciently small domains. 

A resume of the results for linear elliptic 
and parabolic equations is given in the opening 
sections of the respective chapters. Goursat’s 
‘‘Cours d’ Anlyse Mathematique” vol. 3, or Frank- 
von Mises’ ‘‘Differential-und-Integral gleichungen 
der Mathematischen Physik” vol. i can be consulted 
for proofs and other details of these results. 

No account is given of the method of Green’s 
function for the hyperbolic equations because it 
is almost similar to the method for elliptic and 
parabolic equations. Besides, an almost exhaustive 
account of hyperbolic equations is given in Courant 
and Hilbert’s admirable book ‘"Methoden der 
mathematischen Physik” vol. 2, recently published. 



CHAPTER I 


NON-LINEAR PARABOLIC EQUATIONS 


1*1. Convergence of an Infinite Series 


Let p(x) be a function > i in the interval 
o < jt, and let p{x) and ^be continuous and 


uniformly bounded in this interval : 




< M, 


\dp 

\dx 


< M, 


M being an absolute Constant. 

Let the Sturm — LiouvUle equation 

« 

for the boundary conditions 

^(o) = o, 9i(3t) = o (3) 

have the eigenvalues A,„ and the eigenfunctions 
4>ni^\ where « = i, 2, ... 00. We assume that the 
set of functions is orthonormal. The asymp- 
totic expansions of and 4>J>d) are known* to be 


♦Curant— Helbert : “Methoden der Mathematischen 
Physik”, Vol. I, Chap. 5, § ii, p. 291. 
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- «%2 + 0(0. « - j WM 


n {“ dx 




nn cos 


0(i).^=j 




« •'oV/)(>^) 


Evidently -3- does not approach zero whatever 

n may he, so that — is different from zero, and there- 
in 

fore bounded for all n. Consequently, <i>Px) 

and — are also uniformly bounded. 
n dx ^ 


Hence we can write for all x in o < < jr 

« 

(« = I, 2, ... ) 

where h is an absolute constant. 

Finally, let the function p^^,{x, t) be defined 
in the domain o<x<3t, o</<oo for any 

r,r > I, and suppose that/„(x, t), are 

uniformly bounded in the whole domain for all 
r, J ^ I : 
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3 


Pr.,i^y t) 


<N, 


^Pr., 


<N, 


^^Pr., 


<-N-, (6) 


N being also an absolute constant. 
Consider the function 


where », A , r, s can take all positive 

integral values from i to 00, and where dashes are 
used to denote differentiation w. r. to x. 


We shall establish the following result : 
For all Ai, A^, /j , 4 ^ i we have 






( 8 ) 


where is an absolute constant independent of 
r, s, and polynomials in r and s 

respectively. 


To prove this, we integrate (7) by parts, and 
get 



(9) 


2 
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But wc have 

= PPr.,-^{K-<l>kryii-^'l, 


+ PPr, • ^*1 • 


• #:.) 


+ PPr. * ^*1 

••• 

8 


~ Pr.$ " 




+ Pr. * ^hi • 



^'s} 

+ PPr. t ^*1 

••• ••• ^ 

V 

(lo) 

We must now differentiate these expressions. 

Thus 

^[Pr.^^h 



■^}] 

= Pr.,Vl,- 
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5 


+ A. •” *^4 

~ ^ti ••• *'* 

-A. i{#., V--f.,i[^:,-'f:.-. 

~ A.t ^*1 ••• ^r ‘P'li ••• ^I.(^Ai + ^*2+ ••• + 

+ PPr ., ••• ^*<+1 •’• 


■■■ *^*>-)’ 

where means that the term / = / is omitted. 

We see that the number of terms on the right 
of (ii) is 

= r -f- arr + r + r (r — i) = r(r + at + i). (la) 

We differentiate now the second line on 
the right of (lo). On account of the relation 
p,f," -f- p' this term is 
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— Pr,» ^ki"’ 'f‘kr ^h'” ^lj-1 

Pr.t ^*1 ••• ••• (^3) 

On differentiating the first line in (13) we get 

= - A. A, - K *‘1 - ^0 


- ■fi, - A, fi;,,. - '‘:.2 (^.1 - 

^*i+i •■■ A,)| 

~ A,. Ai ••• 'i’k^ i>'ii ••• ••• + \) 

-A,, ■f.,- - *‘.-1 

■■■ ^Ol (*4) 


The last term on the right of (14) is equal to 


^kl ••• P' ^0 


<t>Lr - <!>', 
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The number of terms on the right of (14) is 
therefore 

= r + + 2.s(s — 1) s = s (r is). (16) 

We differentiate now the second line in (i3)j 
and get 

^P' Pr.t ^*1 ••• 

— ^P"Pr,t ‘f’ki'" "Aj, ^P' Pr,, 

•f>kr 'Aji 


— ^P'Pr., ... 4>k^) 

- ^P'Pr.. K 

(^7) 

The last term on the right of (17) is equal to 

+ P ' ^'0 

••• 

Thus the number of terms on the right of (17) is 
j(2 + r + 2r) (19) 
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Finally we differentiate the last line on the right of 
(lo), and get 


^*1 ••• ^li ••• 



— P*P'r.t^ky"’ ••• PPr., ”• 

+ PP't,A\x ••• ••• ••• ^*r) 




The last term on the right of (20) is equal to 


Pr, t^k-y” 




^tj +P''f’'tj) 




Thus the number of terms on the right of (20) is 

-= 2 + r + 2J. (22) 

Therefore the integrand in (9) has altogether 
m terms where on account of (12), (16), (19) and 
(22) we find 

ffi = r(r -{- zs i) + s(r + 2j) + ^(2 + r + 2r) 

+ (2 + r + 2j) = (r + 2r + i)^ + I (23) 

Moreover, each term in the integrand contains 
(a) either , (x, t) or X. . (-v, t) or X! , 0 > 
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t/ 

(b) either p{x) or p%x) or^ 


/2 

or-s- or I 


(c) either or i (y = i, 2, ... r) ; 

{d) either or i (y = i, 2, ... r) ; 

{e) either <}>^^ix) or {x) (y = i, 2, ... r) ; 
(/) either or <f>\. {x) (y = i, 2, ... r) ; 


{&) 


To make the result independent of k^, 
we must divide the expression by 
one and the same term may contain and to- 
gether, therefore it is not sufficient to divide by 
'K^^ only, but by Thus the whole expression 
must be divided by (A,j^ ... Aj^)^''-. Then each term 
in the integrand will be uniformly bounded for all 

», ..., /i, ..., /, ^ I. 


We shall denote the m integrals thus arising by 

(24) 

From (9) we obtain then 
VA„ \a\ . I 


-- C V UiA 


(»0 


From Parseval’s theorem for the Sturm- 
Liouville eigenfunctions we know that the series 
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2 convergent for all / = i, 2, ;5??, so 

n 

that 

2 ^2(r+.) (26) 

n 

From (25) we obtain on summing over n and 
applying Schwarz’s inequality as well as the inequa- 
lity zab ^ : 


< 2 . ^27) 

where we have used the result that 


Writing 


^ \ ^ " w2 

n 


6 • 


N^M\ 


and 


^1 {r)=^r+ I, 

% (j) = 2J + I, 

the truth of (8) becomes obvious. 


J, 


(28) 

(29) 


i‘2. An Infinite System of Non-linear 
Integral Equations. 


Let the numbers «, k^, /i, •••, /,> r, s 
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take as before all positive integral values from i 
to 00, and let /„(x), gjx, j), hJ^j, ky, ly, 

I, r,s) be sequences of functions, denoted 
simply by /„, g^, defined for all values of x,j : 

o < < 00 ; o < x" < 00, ... (i) 

and satisfying the relations 

2 I fn{^) I ^ ^ fot all X in (i) ; (2) 

1 j) 4 y fot all in (i) and all « > i; (3) 
Mo 

'L\h^\ < A q^if) q^s) for all y in (i) and 

n ' 

all and /’r > i. (4) 

A, B, b, c are absolute constants, and qy{r)y 
q^(s) denote polynomials in r and s respectively. 

We consider the infinite system of non-linear 
integral equations : 

»n('V)=/»M+r^n(^'-.j)^ ^ K»ky{.y)- 

r, s k^, kj. 

4, 

% 0 ') -"jiO') •••"(.( (5) 

(« = I, 2, ... go) 

We shall prove that one and only solution of 
the system (5) exists for aU values of x, provided 
only that e is less than a certain constant. 
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We shall solve the system (5) by the method 
of successive approximations, and write for this 
purpose 

(x) =/„(x), (6) 

and for all /!!? > i 

(^) j) 2: K - 

k, 

(j) (j /) dj (7) 

We must prove that these approximations are 
convergent, i.e., that lim„_^^ (x) exists for all 

« > I. This will be the case if the doubly infinite 

00 

series ^ ^ | (x) — (x) | is uniformly con- 

w=0 n 

vergent for all x. 

We begin by proving that the series 

( 8 ) 

n 

is uniformly convergent for all x and all m. 

Summing (7) formally over «, and taking ac- 
count of (2), (3) and (4), we obtain for all x 

2 I (x) I < ^- + d 2 qx{r){b Max 2 | I'X 

n r n 

2 Max 2 [ (^) | }*, (9) 

8 n 

where we have written = AB. 
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Since qi(f) ana q^{s) are polynomials, two constants 

d and e* can be found such that 

^i(r) < d^ for aU r ^ i, fot all J > i (10) 

From (6) we get for all x > o 

= ... (IX) 


Substituting (10) and (ii) in (9) £ot m = i, we get 
2 ^ \'L(hcdy. 

n y r s ) 

Let f< Min then 

2 + t^'k- 


If c is taken to be less than 

z , bed bee 

Wd'e+ bd + be ’ i bed V-We ^ 

and we obtain 

2 ( <'’(x) \<e + e= ze. (12) 

n 


Writing m = z in (9), and substituting (12), we get 

2 e + a 1 2 (zbedy 2 (ibee)" I . 

n V r # y 


*e must not be confused with the base of the exponential 
functions. 
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Assuming that c < Min 
_ ,„w . ^ . ibcd zbce 

Supposing that 

^ ^ zb'^de ^ bd be' 

0 h/'/i ‘y h/'fi 

which makes a , , — , < we obtain 

1 — zbcd I — zbce 

2 1 I < ^ 

n 

Thus we find generally for all m \ and all x ^ o 

2 I u^T\^') \<zc, (14) 

n 

establishing the uniform convergence of the series 

( 8 ). 

We shall prove now that the double series 

2 2 (15) 

?n==0 n 

is uniformly convergent for all x ^ o. 

In fact we have from (7) 

= I j) ^ ^ K 

•^0 rs k\ fcy, 

ll Is 

X (jy) - «i:> (j) - ^r’Xj'Kr" (j) 

- (j)] 4 y 
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= ^ 2 : h{u[”^\y) ... 

Jq ky> V 









(,,<»> <j) - 

^(m— l) 

O’)) 

+ 


(j) 

«(«) 

(j)-- 

J) «10)- 


O)) 

+ 

£/(m— 1) 

(j) 


(j) 

"iir* (j) ■•• (j) («i”\ 

>0) 






- 

^a-2 

(j)) 

+ 








+ 

^(»n— 1) 

(J) 


■” 0) 

... «<;-» (:>')«' 

(J) 



- «{r‘’ w)}^- 

The uniform convergence of ^ | (x) — (x) ( 

n 

is seen therefore to depend on that of 'L ( 

n 

— (x) |, so that on account of ( 14 ) we obtain 

for all X 

2 (^) ~ (•’^'■) I ^ + -t) 

n r.« 

% W X Max 2 1 (x) — (x) I 

n 

< a 2 (r + i) (y-cy-' qy{r) X 2 (j + i) ^ 

r a 

{zc)’ (j) X Max 2 ( »1"*' (x) — (x) ( 

(16) 
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f 

But we can find two numbers and such that 

(r + i) ^i(r) < for all r > o) , 

(s + i) ^2^) ^ fo^ all j ^ o J '■ ' 

The relation (i6) then becomes 
2 1 (x) — (at) \<ab^ {zbcd^'-. 2 {zbce-^‘ 


n r=() 8 

X Max 2 1 (>r) — I (i8) 

n 


If we assume that ^■< Min(-^ , — ), we get 


2 I (x) — u {x) I < ab 

n 


zbce-^ 


(i — zbcd^ (i — zbce-^ 
X Max 2 ( {x) — (a-) I (19) 


Finally we assume that 

^ < XabeiVd:,-^~i^ ’ 

which gives 

zab"^ ce-i . , 

^ ” (i — zbcd^(i — zbce^) ^ 

Repeating the reduction formula (19) m times, 
we obtain 

2 I (x) — ul”'^ (x) [ 

fl 

< Y’" Max 2 I »n^ (x) — uy (x) [ (zz) 
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Summing over m from o to 00, we get for all x > o 

Z Z I (x) - wl”) (x) I 

m*>0 n 

< Z Y”* Max Z I (x) — (x) I (23) 

m n 


QO 


Since Y < I, Z Y”* is convergent. Moreover, since 

w==0 ^ 

Z I (x) ( < on account of (ii), we find for 


n 

all X > o 


2 1 (x) — (x) I < Z 2 

n r 8 

bed bee 
^ ^ I — bed ' j — bee’ 


(^4) 


since both < i and bee < i. 


From (23) we deduce therefore that the doubly 
infinite series (15) is uniformly convergent for all 
X > o. 

From this we conclude that all the limit func- 
tions 

^»i^) ai”’ (^), (« = 1 . 2. 3 » • • •). (25) 

tn -> 00 

exist and are continuous functions of x. 


Moreover for all x we obtain 
2 1 »„ (^) I = 2 Urn I (x) I < ze. 

n n m-^00 


(26) 
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From (7) we get, on passing to the limit 
»n (^) = In W + r (^, j) 2 2 

•'O r.a Afj. 

Al. . .Is 

X % (j) . . . (j) (j) . . . {j)dy, (27) 

showing that the function {x) obtained from 
(25) is a solution of the integral equation (5). 

We have to show now the uniqueness of the 
solution, under the condition that 2 ( (•v)| 

n 

should be uniformly convergent. 

If possible suppose that the system (5) has 
another solution (x), {n = i, 2, ...) such that 

2 [ (^) 1 < (28) 

Then 

Vn M = /„ (^) + j" gn (^> J) 2 2 h„ 

•'O r.« Acj, 

l\, . Iff 

X % ( j) ( j) (j) . . . Vi^ {y)dy (29) 

Then we get from (7) and (29) 

M M = 1 .gn j) 2 2 

‘1 ‘s 

(j) ••• (j) - w 
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X {% (j) [\ (j) - O)] 

+ Vki (j) (j) (j) - 

+ 

+ o) o) u) 

KiO)- 4 r‘’( 7 )]} (30) 


On account of (14) and (28) we obtain therefore 
^ 1 (x) 1 < >: (r -b j) (r) ^2 (j) 

n r,<? 

(2c)r+s-l X Max I V (.v) — (a") I 


zab^re^ 

(l — 2/Wj)(l — 2/Wj) 


Max 2 1 (x) — (x)(, 

n 


as found in (19). Repeating this inequality ;;; times, 
we obtain for all x > o 


Z [ (x) - /4’”> (x) [ < Y”’. 

X Max r, 1 (x) - (x) ( (31) 

n 

where y is defined in (21). But 



>' %U) \ O') \ U) \ (y) 4y^ 


so that we get on account of (28) for all x ^ o. 




O) — »n ^ (x)\<a:^ b'-+‘ (r) q^ (r)( 2 <r)’-+' 


<a 


zbed ibce 
— zbed I — zbee ’ 


( 3 ^) 
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since (r) < and (s) < 

Therefore, since y < i, we get from (31) 


lim 2 [ {x) — (x) I = o (33) 

m ->• 00 

Consequently, for all x > o 

vM= ( 34 ) 

m -> QO 

(« = I, 2, 3 , • • • • 00), 


showing that the two solutions are identical. 

1.3. Solution of the first boundary problem. 

Let pix) be a function > i in the interval 
o < X < A, and let p{x) and dpjdx be continuous 
and uniformly bounded in the whole interval. 


We consider the non-linear differential equation 


« 


d 

6 x 

in the domain 
o < X 


3t, O < / < 00. 


(2) 


We shall prove that there is one and only one 
regular solution of (i) which satisfies the conditions 

u (o, t) = u ( a , t) ~ o for all / in (2), (3) 


and 


u (x, o) — /(x) for all x in (2). 


( 4 ) 
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We assume that the given boundary function 
f(x) can be expanded in a series of sturm-Liouville 
eigenfunctions, i.e., the eigenfunctions of the 
equation 

+ I (5) 

q)(o) — <p(it) — o. j 

Thus 

Ax) = ^ ( 6 ) 

n 

where the coefficients are such that the series 
S is absolutely convergent, being the 

n 

eigenvalue of (5) corresponding to tp^. Hence 

2 k J (7) 

n 

where c is an absolute constant. 

The question of the expansion of a function in 
a series of the eigenfunctions 9 „(a'') has been 
considered by Lichtenstein.* 

The function P(,v, is supposed to be 

analytic in the variables u, , so that it can be 

* L. Lichtenstein : “ Zur Analysis der unendich vielen 

Variabeln”. Palermo Renclkonti 38 (1914), 113-166. 
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expanded in a power series 

^ ^ dx) "" ( dx) • • • • 

The coefficients and their first two 

^ d ‘^p 

derivatives and * are continuous and uni- 

ox dx^ 

formly bounded in the domain (2). For values 
of I » [ < I, < I, the series (8) will then converge 

absolutely and uniformly in the whole domain (2). 

The solution of the differential equation (i) 
will be determined as a series of the eigenfunctions 

<p ; 

«(a-, /) = X 1 ’„ (/) <p„(x) (9) 

n 

The conditions (4) and (6) require that the 
coefficient p„ (/) should have the initial value : 

•••) 

Now since ‘P„(.v) is a complete set of ortho- 
normal functions, we have 

dx) (“) 

where 

^„(r, J ; /) = S ^n(^> A> •••> ^') 


k] k^ 


X %(0 . • • %(f) ■ • • (1 i) 
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and = a^{t, k^, . . k^, /„ r, s) 

Jq i > 1 8 

(13) 

dashes denoting derivatives with respect to x. 

Assuming for the present, what will be proved 
in the sequal, that the two series 

dv 

2 and ^ " tp„(^) (14) 

n n 

are absolutely and uniformly convergent in the 
domain (2), we get on substituting (9) and (ii) in 

(1): 

- 2 ^ - 2 -- 2 ^ Zn(r^ 

n n n r ,s 

(13) 

For each « > i we get therefore 

-h K i’nQ) = - 2 -f ; 0, (16) 

of which the solution satisfying the initial condition 
(10), viz. F„ (o) = c„, is 

i'nW — ( e-^n<f-y^ S ^„(r, s; y) dj 

*'0 S 

= 2 2 a,v^{y) ... 

0 r, s 1(1 ^ .... ky ^ 

Ih •••• ^8 

% ( J>ji( j) • • • \{j) dy, . . (17) 

« = I, 2, ... 00. 



24 


PARABOLIC EQUATIONS I 


This is an infinite system of non-linear integral 
equations for the determination of the coefficients 
vj^t) in the solution (9). 

We write for all « > i and all / > o : 


» n ( f ) = (18) 

/n (0 = (19) 

(20) 

•••5 ^r> *•*5 ^s’ ‘^) 

(21) 

('*! ••• h,h. ••• 


Then multiplying the equation (17) by we 
obtain : 

»«(0 "= /n(0 -I- (' SJJ, y)^^ ^■’n%(y) ••• 

•'O r. s rci fcy. 

ii is 

\(y) »i,(y) • • • »iSy) 4y^ 

(« = 1, 2, ... 00). 

This is the same system of integral equations 
as considered in § i * 2, except that we have now 
written t for x. 

In order to be able to use the results of that 
section, we must prove that the functions /„, .g„, 
as defined in (19), (20) and (21) satisfy the relations 
(2), (3) and (4) of § I • 2. 
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» 

Now, for all / ^ o we have, on account of (7) : 

S I /»(0 1 = 2 ( 1 < 2 1 itK f < ( 25 ) 

n n n 

Further 

^ gnit^j) <b = — K\ e-^ntVeKy dy 

Jo Jo Jo 


so that for all « > i and all / ^ o, we have 

I I < I ( 24 ) 

Jo 

Also the function />„ satisfies the relation (4) 
§1.2, as shown in § i ’ i. 


Therefore, as proved in § i • 2, the sytem of 
non-linear integral equations (22) has one and 
only one solution {t\ given for all t by 
«« if) "i”’ (4 (« -- I, 2, ... 00), ... (25) 

where (/) is the successive approximation 
defined by 

and for m ^ x 

(f) = fn (^) + r ^ ^ 

•'O r,8 k\, ky, 

h h 

»ir“ ( a) • • • o) ( j) • • • 4y- w 
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It was further proved that 2 [ if) [ is uni> 

n 

formly convergent for all / > o. 

If we now take the coefficients (/) given by 

K {t) = ~/2 «„ (^), (« -- I, 2, ... co), (27) 

then (/) satisfies the integral equation (17) and 
consequently the differential equation (16). Further, 
on substituting the value of (/) given by (27) in 
2 2 [ .f ; /) [ , this series is easily verified 

n r.s 

to be uniformly convergent for all /, on account 
of the fact that 2 | {t) | --- 2 | {i) | has the 

n n 

same property. Also, since the functions (.v) 
are uniformly bounded for all a’ in (2) and for all 
« > 1, the series 2 (/) q)„ (a) and 2 2 

n n IS, 

dv 

f -f ; and consequently 2 (x) 

are absolutely and uniformly convergent in the 
domain (2). 

Thus 

u (a, /) = 2 (/) 'P„ (a) (28) 

n 

is the one and only solution of the differential 
equation (i) satisfying the boundary conditions (5) 
and (4). 
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1.4. An Equation of the higher Order. 

In §1.3 we considered an equation of the 
second order. In this section we shall extend the 
proof to an equation of 4th, 6th, or any higher 
even order. 


For convenience, we take the case when pix) 
= I, so that the sturm-Liouville eigenfunctions 
<P„(,v) degenerate into the sine functions sin nx. 


Let m be any positive integer > 2, and let the 
equation 




(~ i)’” 


hti 

a/ 


P (x, t ; u. 


\ 


(i) 


be considered in the domain o</<oo 

and let it be required to determine a regular solu- 
tion, satisfying the boundary conditions : 

a(o, /) ^ : u (n, /) =- o for all / ; (2) 

u(x, o) /(x) for all .v in o < x < re. (3) 
The function f(x) can be expanded in a series 
/(x) = 2 sin »x, 

n 

such that 2 [ is convergent* : (4) 


^2m f 

* For this it is sufficient that exists, and can be re- 

presented as an indefinite integral, cf. E. W. Hobson : 
Theory of function of a real variable^ vol. 2, p. 516. 
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^n^”'\e„\ = c. ( 5 ) 

n 

The function P satisfies the same conditions 
as in §1.3. 

Thus P (x, / ; ) = 2 p,. . {x, t)if (— )‘ (6) 

The solution of the equation (i) is then deter- 
mined as 

u{x, /) = 2 i'n (/) sin «x, (7) 

n 

so that on account of (3) and (4) we get 

?’n(o) = ^« («= I. 2 , ...) (8) 

Then 

Pr. . (•^> ^ ; 0 sin «.Y, (9) 

where 

Zn{r, -f ; 0 = 2 «n %(^) • • • ^(,(0 • • • 

k\ ky. 

1 1 Is 

/j, /g, r, J*,) 

I TT 

“ 1 Pr • sin kyX ... sin cos l^x ... 

Jo 

cos /.x sin nx dx. (i i) 

Substituting (6), (7) and (9) in (i), we obtain 

( — i)”* 2 iP”' vjf) sin + ( — i)”* 2 ~P sin nx 

n 71 

= 22 ^„(r, J ; /) sin «x, ( 12 ) 

n r, 8 
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where all tlie series in (12) have been assumed 
to be absolutely and uniformly convergent. This 
will be proved later. 

From (12) we get for each « > i : 

dv 

vjt) = ( - I)" ^ (r, r; /), ... ( 13 ) 

of which the solution satisfying (8) is 


Vn{t) ‘ + (— i)” (* 2 K.n{r, s ; j) dy 

•'0 r. 8 

+(— i)”* 

*^0 r, 8 fc] ky. 

^8 

X J) • ■ • \iy) (14) 

We set for all / ^ o and all ^ I : 


»„(/)- - /-„(/) ; ( 15 ) 

/« (0 ; (16) 

(A j') (- !)”■ ; (yi) 

K=K 0 ; ^ 1 , ■••> •••. 4 ; '■> -f) 


-‘( k ,.: rtX -- i3 ^ ; 

From (14) we obtain, on multiplying both sides by 


„%m 


»«(0 


0 >“.8 ’^1 fcr 

h* •••• ^8 


S ( J >') • • • ( J ) "(1 O ') -\ U ) 4y > 9 ) 

(n = I, 2 , ... 00 ). 
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It can be proved as before that fl, g„, satisfy 
inequalities (2) (3) and (4) of §1.2. Hence the 
infinite system of non-linear integral equations can 
be solved uniquely by the method of successive 
approximations. From (/) we can then deter- 
mine if) according to the relations 

Vn if) = (/) (« = I, 2, ...) (20) 

Since 2 | f) j is uniformly convergent for all 

n 

/, it can be easily verified that all the series in (12) 
are absolutely and uniformly convergent. 

Hence 

«(.v, t) - -- 2 f) sin nx 

n 

is the required solution of (i) satisfying the boun- 
dary conditions (2) and (3), and this solution is 
unique. 


1.5. An Elliptic-Parabolic Equation. 


We consider the differential equation 

bx^ aj2 ~ d/ “ " ’ 

and determine a solution uix, y ; t) which is regular 
in the domain R defined by 

o < X < Jt, o < j/ < 31 , o < / < 00, 


(2) 
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and which satisfies the boundary conditions : 

»(o, y, t) u(K,y ; /) = o for all y and / in R ; K . 
u{x, o ; /) = u(x\ jc; /) — o for all x and / in R ; J 
and 

u{x, y ; o) =/(-v, )') for all a- andj in R (4) 


We assume that can be expanded in the 

double Fourier scries 

f(x,j) sin ;;/.v sin (5) 

m,n 

and that the series fi- n^) ( [ is convergent. 

1)1, n 

For the solution of (i) we write 
u{x, y ; /) v,„_„ (/) sin mx smny (6) 

m,ii 

and get 

u\x, y ; /) -= (/) sin mx sin ny, (7) 

VI, n 

whe re 


(0 “ ^ /"« v) (/) (t), 

K\fJ.U 

with 

'y TT 


fm ( 5 ^. v) 


K . 


sin x-v sin lx sin mx dx. 


0 


fn (i‘> ^ ! I ^ 

Jl Jq 

From (i), (6) and (7) we get then 

— i. + «2) (/) sin mx sin ny 


( 8 ) 

(9) 



32 


PARABOLIC EQUATIONS I 


— 2 sin ny 

= 2 sin wx sin ; (lo) 

therefore for all \ 

+ (w2 + «2) ^ ^ (l l) 


We have assumed that all the series in (lo) 
are absolutely and uniformly convergent in the 
domain R. As usual, this can be proved later. 

From (4) and (5) we see that the solution (6) 
would have to satisfy the initial condition 

V„,n (o) = Cr,^.n (W, « 1 , 2, . . . ) (l z) 

The solution of (ii) which satisfies (12) is 
given by ; 

Vn (0 - I' 

•^0 


— *'ni,n ^ \ ^ ^ J /7 

•^0 h\fiv 

(S)t\ , {S)ds (13) 

We write for all n'^ \ : 


»'«.n W = d- «^) ^m.« W, d- «^) 
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and get from (13) 

•'O 

X y n, (s\nj (s)ds flO 

This is a doubly infinite system of non-linear 
integral equations for the determination of the 
functions ^ (/) which in their turn determine the 
Fourier coefficients „ (/). We shall solve this 
system by the method of successive approximations, 
and for this purpose we set 

K]n Q) 4.n (l6) 

and for all r > i 

W = ^'m.n + tfi) j f-<«2+n2)n-.) 


X 




We have to prove the convergence of this 
approximation. We shall first show that the 
doubly infinite series 


„ X + !>”) (l* + V*) 

is uniformly convergent for all x?.(iv. 


(18) 
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From ( 9 ) we have, on writing A 


d2 52 

dx^ ^ dj/2 ‘ 


/. (X. »/. (C, v) = - (' j' (sin 

sin \i y sin "kx sin vj) x A (sin sin nj) dx dy 

X A (sin XAT sin \iy sin kx sin vj/) dxdj 


= - I' C V 

{— (x^ + k^) sin XA- siagj' X sin kx sin vy 
+ 2xX cos x.v cos kx sin \iy sin v j ■— (jx^ [- v-) sin 
XX sin ny sin kx sin vy 4 - 2[iv sin xx sin kx cos g;/ 
cos vy}dx dy. 


Therefore 

|/m(^> ^)/n (l^s ■'') I ^ 4 i \ f (yi k') f (ll v'i I 

(X2 + |a2)A2 v2) ^ («?- - I- «2) V I 

+ ^)/„(^, V) + |/„,(X, X)^„(g, V) I I , (19) 

where 

2, 

gmi^i “ 1 cos XX cos kx sin mx dx, 
n Jq 

^„(g, v) = — [ cos (ix cos vx sin nx dx. 

X J n 


^ (20) 
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From Parsevars theorem we know that for all 
Jt, Jl, I*, V : 

S/2(5 <,X)< 2, 

n n 

2 X) < 2 , 2 ^^|X, V) < 2 . 

n n 

Therefore 

2 v) < 4, 'Lgl{v., X)/^(n, v) < 4, 

m, n ?», n 

2 fUy-, v) < 4. 

m, n 

We get from (19) on summing, squaring and 
using Schwaz’s inequality as well as the inequality 
zab < : 

I 1 /m(x, I f < I N 1 

+ v2)) ^ 4' „2)2) 


2 

m, n 


^)/„(|A, v) I -f [gj%, l)fSV; v) 

+ fl 


+ >fn(«. ^)/n(^ V) -f /|.(K, X)^^(h,, V)| 

^ V 3J< 


But the doubly infinite series 2 > - - g f — 2^ is 

m, n \P^ I ^ / 


4 
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convergent, so that 

V)h 2 


and with it the series (18) is uniformly convergent 
for all X, )i, V : 


^ !/».(«> ^)/ n (t*» V ) I ^ 

«^«(»t=*+|X^)(X2+v2) 


(21) 


where <7 is an absolute constant. 

We write, since 2 + «^) n is assumed 

m. n 

convergent, 

2 I == 2 {nP- + tP) (22) 

m, n m» n 


We have also for all / > o and all «?, « > i: 


11 ' 

Jo 




c— (w^d- 


n2)i r 

Jo 




X j 


I (/ W2 + « 2 ) < (^2 + „ 2 ) • 


(23) 


Thus from (21) and (25) we get for all / > o 
and all x, A,, n, v > i : 


w, n Jo 


^)/«((*. V) I 
(x^ + + P) 


< 2 + fP) • 

w, n 


I ^)/n(^» V) I 

(/»* + ■ (x^ -j- + P) 


<a. 


(24) 
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Then we can show as in § i . 2, that the approxi- 
mations (17) are convergent, and that the system of 
non-linear integral equations (13) has one and only 
one solution. The only alteration in the proce- 
dure would be to change the simple summation 
over « to a double summation over m and n. Thus 
it is seen that the non-linear elliptic-parabolic 
equation (i) has a unique solution satisfying the 
boundary conditions (3) and (4). 

1.6. The mixed boundary problem. 

In the previous sections we considered only 
the first boundary problem in whch u had given 
values at both extremes of the x-interval. The 
method is, however, capable of dealing with other 
kinds of boundary problems. We shall consider 
here the third boundary problem, the so called 
mixed problem, in which the value of « is given 

at one end of the x-interval, and that of ~ is given 
at the other end. 


Thus we try to determine a solution »(x, /) of 






r 


(0 


in the domain 

o < X < «, o < / < 00, (2) 
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such that u(x, t) satisfies the mixed conditions 


«(o, t) — o, ^ (jt, /) = o for all / ^ o, 
and 

u (x, o) =f(x) for all >r in o < X < 3t. 


(3) 

(4) 


This problem can be brought to depend upon 
the boundary problem we have hitherto considered 
by continuing the functions u(x, t), {x, /) and f(x) 

in the interval n < < 2 Jt according to the method 

of images. Thus we take for all / > o : 

U (231 — X, /) = u {x, f), j 

(23t - x,t)=p^{x,t), ^ ( 5 ) 

/(23t - x) =f{x), J 
where varies in the interval o < x < :t. 


Now we define functions (x, t), p'^ (x, t) 
and /i (x) in the interval o < x- < 23t, for all / > o 
as follows : 


^1 (•^> 0 = /) in o < X < 7t, 

— u(zn — X, /) in :t < X < 2Jt; 

Pr (^> f) = Pr{^^ /) in O < X < 3t, 

= pjiz'n — X, /) in 31 X < 23t; 

/i (•^) =/(^) in O ^ X ^ 31, 

= /(2Jt — x) in 3t < X < 23t. 


( 6 ) 

(7) 


(8) 
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Then we see*that % (x, /) is a regular solution of 


d^u du / A / N 

= W 

in the domain 

O < X < 2rt, O < / < CXD, (lo) 


which satisfies the boundary conditions : 

//j (o, /) = « (o, /) = o, I 

/) = a(2Jt — 23 t, /) = a(o, r) = o, j ' ' 

and 

u-^(x, o) = »(x, o) — /(at) in o < at < Jt, 

Ui(x, 6) — u{zn — X, o) =f{zn — x) in ji<x<2a, 
so that 

Uy(x, o) = fi{x) in o < x < 2re (12) 

But the method of § i . 3 is applicable to an 
interval of any length, and therefore the solution 
«i(x, /) of the equation (9) is unique. 

We define another function U 2 {x, t) in the 
domain (10) by the relation 

u^{x, t) = u^(zn — X, t). ( 15 ) 

Then we get for all .v and t in the domain (10) : 

92^2 {X, t) 9^2 (X, t) ___ dhiy (2Jt — X, t) 

9x^ 9/ dx'^ 

9aj (2Jt — X, /) 

~ 9/ 
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= 2 pri^n — X, — Xy t) 

r 

= 'Z p'riXy t)ul{x, t), ... (14) 

r 

showing that u^ix, /) satisfies the diflferential equa- 
tion (9). Further 

«2(2n, /) = «i(o, /) = o ; J 

and for all in o < X < 2 Jt 

u^ix, o) = Ui(zn — X, 6 ) - x) =^fi(x). 

(16) 

Equations (15) and (16) signify that U2{xyt) 
satisfies the same boundary conditions as «i(x, t). 
It must therefore be identical with u-y{x, t) : 


i.e., on account of (13), 

«i(2Ji — Xyt) = «i(x, t) (17) 

Hence 

= ° (■») 


Thus, in the domain (2) the function «i(x, t) 
satisfies the differential equation (i) and the boun- 
dary conditions (3) and (4). It is therefore the 
required solution. 
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1.7. Non-rectangular Domain. 

So far we have considered only rectangular 
domains 

T, (T-^cx)) (I) 

Now we take a general domain D bounded by 
the x-axis between the points o and n, and the 
two curves Q and Cg with the equations 

(Qx = i,(/) (2) 

where 

hi (o) = o, h^ (o) = n, ^2 (/) > hi (/) for all /. (3) 

The problem is to determine the solution 
«(x, /) of the differential equation 

(4) 

which is regular in D, and which satisfies the 
boundary conditions : 

« = o on Q and on Cg, ) / \ 

« = J{x) for f = o. 1 

We shall prove that on transforming to new 
variables, the problem can be reduced to the 

determination of the solution of 

d^u ^ 


(6) 
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in a rectangular domain of the typef(i), and for 
analogous boundary conditions. Our integration 
method as developed in the previous sections 
would then be applicable to it. 

For this purpose we first make the transforma- 


tion 

X = (/) x’ g {t)jn, / = /, (7) 

where 

= ( 8 ) 


To the domain D of the (x, t) plane would then 
correspond the rectangular domain 

O < x' < Jt, o < (9) 

of the (x'j /) plane. The differential equation (4) is 
then transformed to the equation 

ft = 2 f) - (10) 


If we transform further 
' ‘ dt 






and write 


t/{x\ /') = »(x, /), g'it') == g{t), 1 

pX^\n^pXx,t),f\x')^f{x\] 


(11) 

(12) 
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then we see that u- satisfies the difFetential equation 


fix'* d/' 




(15) 


and the boundary conditions become : 
a' (o, /') — u'(n, /') = o for all / > o, 1 
t/(x\ o) for all x' in o < x' < it. j 


(14) 


If we write 

o = -'4 o. - (■!) 

and drop suffixes throughout, the equation (13) is 
seen to be the same as (6). Further, it is obvious 
from (ii) that 00 with /. 

The proposed reduction has thus been obtained. 



CHAPTER II 


NON-LINEAR HYPERBOLIC EQUATIONS 

2*1. Solution of a system of non-linear 
integral equations in the restricted domain. 

We consider now an infinite system of non- 


linear integral equations 

»'» (^) = In W + r gn j) 2 

•"o h. *2 hr 

"'*1 (j) »'*2 O) - {J)dy> (l) 

(« = I, 2, ... 00 ), where /„ {x\g„ {x,y) and (j, 
ki, ..., r) are sequences defined for the values 
of X, j in the domain 

o < A" < T, o < J/ < X < T, (z) 

and satisfying the relations 

2 (/„ (^) [ < for all X in (2) ; (3) 

” Un (^,j) I < for all X, J in (2) 

and all « > I ; (4) 

2 1 1 < ^ h' pif), for all y in (i) 

n 

and all k^y > i (5) 
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Ay By c are again absolute constants, and 
pif) denotes a polynomial in r. It will be seen 
that the functions and satisfy inequalities 
similar to those in § i ’2, while (•Vjj) itself is 
uniformly bounded now, whereas its integral was 
assumed to be uniformly bounded in § i ‘2. This 
fact brings about the fundamental difference that 
whereas the solution of the system (5) § i ’2 existed 
for all values of x however large, the solution of the 
system (i) will exist only for a restricted domain x. 

We transform the equation (i) by writing 


K M (^) -/» (^). (« = I. •••)> (6) 

so that (i) becomes 

(x) = {^,j) 2 b„ (/j (j) + (j)) ... 

Jo ky..,»kr 

U'krO) + \{jy)4y> ( 7 ) 

(« = I, 2, ... 00). 


We use again the method of successive ap- 
proximations to solve the system (7), and write 

(x) = o, and for all m ^ i : 

(x) = r g„ (Xyj) 2 b„ (/j (j) + - 

Jo fcl ^ 

(fkr (j) + U)) 4y- ( 8 ) 
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We have to show now that these approxima- 
tions are convergent. For this purpose we prove 
first that the series 2 [ (a-) [ is uniformly con- 

ft 

vergent for all x in a certain interval, and for all 

m i. 

Let 

— Max 2 ( (x) 1 for all x in (2) (9) 

n- 

{f» = O, I, 2, ... 00). 

Then from (8) we find on account of (3), (4) 


and (5) : 

< ABT h^p{r) {c + (10) 

We define another sequence {R„} as follows : 
Rq = 

2.^-^ a T b'-p{r) {c'- + R;,_i), (i i) 

{m = I, 2 , ... co), 


where we have written AB — a. Now if a and 
jS are any positive real numbers, and if r is any 
integer, we have the inequality (a + 

(a'' -j- yS’’). Thus from (10) and (ii) we see that 
for all O’. 

< Rm, (12) 

so that the sequence is certainly bounded if the 
sequence R„, is bounded. 
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Let q be some fixed positive integer, and let 


for all ^ 

(13) 

Then from (ii) it follows that for m ^ q, 

= 2’-" ^ T b'-p{r) (R;;. - R';,_p 

= /,. T Rm-i + fi-i-i + ••• 

+ Rrix)(I^n.-R.-l), (14) 

where we have written 

P= ab^ p{r). ( 15 ) 

From (14) we get for m ^q\ 

( Rm+i — Rm I < ^ 4 7 ’ 1 R,„ — R„-1 (. ... (16) 


But since 

R,,i - R, + (R2 - Ri) + (R3 - R2) + - 

+ (R«+i “ RJ* 

we have from (16) : 

1 R,.! I < i T.’ (x + (rp T) + K Tr + ... K T)«} 

I Tr , , 

(’ 7 ) 

provided rl^T < 1. 

Now if we select T so that 

- (■*) 
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Then 

R,+i < z/^T^ < I. ( 19 ) 

We find therefore that if R, < i then R,+i < i. 
Since Rj = /, Tf’’ < i, it follows by induction that 
R„ < I, («? = o, I, 2, ... 00) (20) 

From (9) and (ii) we conclude that 

2 1 «« (x) 1 < i, (21) 

n 

for all > I, and all x which lie in the restricted 
domain given by (18). 

We go on to prove now that the double series 
2 2 (x) -»<,”•) (x) I (22) 

?n=50 n 

is also uniformly convergent for all x in o < x < T. 


From (8) we have 

^^M+l) _ ffim) 

= \lgn (^,j) 2 {(/^ + <>)••• (4 + <') 

- (Ai + 


= gn (^>j) 2 2 ...fui-JkUl ••• 

•'0 ...» fej. j“l 




,Am—l) 


)} 


+ 2/^ <-“} 

© 

+ 2/^ ...4_3{<4 - <-2' 

© 
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+ 

(rll) 

+ {<> 

(^5) 

Using the inequalities that 2 [/„ (x) | < c and 

n 

2 1 wi"*’ (>:) I < I for and all at in (2) all ^ ^ i, we 

n 

get from (23) : 

2 I «<”*+“ (x) - «<”•> (x) I < BAb<^ p(r)T { (^) r-i 
+ zQ cr-^ + + ... +(r-l)(^_^J^ 

+ r I Max 2 | wj,”’ (x) — (x) [, 

< „ TMax 2 I (x) — (x) |, (24) 

n 

where 

— + '^ j(^ 5 ) 

We see therefore that the series 2 [ (x) 

n 

— (x) [ is uniformly convergent if the series 

2 [ (x) — I same property. 

n 

Repeating the inequality (25) m tirnes we obtain 
2 1 <“+i’(x) - u^'^'Xx) I < (a,.„ T)™ Max. 2 1 »i“(x) 

(^ 6 ) 
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We take the domain T so small that 



(27) 


where y is any number < i. Then we get from 
(26): 

2 2 I ( < 2 Y”* < (28) 

m=0 n w=0 I Y 


showing that the double series (22) is uniformly 
convergent, as required. 

It follows that all the limit-functions 
u„(x) = lim„ ^ <”‘>(^), (« = I, 2, . . .) (29) 

exist and are continuous. Further, we get on 

account of (21) and (29) : 

(30) 

n 

from (8) we get, on passing to the limit /w -> co, 

y) 2 K{Ky) + u^,{y )) ... 

{fkriy) + »kriy))‘b ( 31 ) 

(« = 1, 2, ... 00). 

Thus ujix), as found from (29), satisfies the 
integral equation (7), therefore wj^x) — «„(x) -j- fjix) 
will satisfy the equation (i). 
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Thus the? existence of a solution of (i) 

is established for sufficiently small T. We shall 
now prove that this is the only solution such that 
2 [ I is uniformly convergent. Obviously, 

n 

it is sufficient for this purpose to show that the 
system (7) has no other solution {x) (n =- I, 2 , ...), 
different from «„(x) as found in (29), which is such 
that 2 [ I converges uniformly. 

n 

If possible, let uj^x) be another solution of 
(7), so that 

K\>^) = [ J) ^ K{ fkyij) + 


( fkriy) + <0')) dy ( 32 ) 

Also let 2 ( uj^x) I < (55) 

n 


for all jv in o < < T, (T being a constant. 

Then from (32) and (8) we find 

=--- f gni^, y) hJ\ ( Ai + K) 

*'0 ky,.,*k,. ' 

( /fc2 + "* 2 ) • • • ( /*r + "fcr) ~ ( /*1 + ) 

( /fe + ) • • • ( A. + (54) 


5 
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From (34) we find as before, for all x in o < a: < T, 
2 I »Jx) — 1 < ^ Max 2 [ Un(x) 

— (35) 

where 

Kfy e,(') = a b^pir) I + (^) Q + i) 

+ Q^-«(^-^ + ^+i)+ ...4 

+ f ... -f ;r4- l) + Qr ^_|_ j)| 

(36) 

Repeating the inequality (55) m times, we obtain 
2 I «„(x) — u^^\x) ( < {^(r, c, i)T\”' X 

Max 2 [»„(x) — u^'^Hx) 1 < c {^(r, c, c)T}”'. (37) 


We assume T to be so small that 

_ JL_ 

p{r, ’ 

where y' is any number < i. Then 
lim 2 I f*Jix) ~~ f^lrKx) j < b lim y'”' 

m-> CO n m-> 00 


” O 
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(39) 


Hence we find that for all in o < x < T : 

= = (40) 
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= ly z, ... ho)y provided only that T satisfies the 
inequalities (i8), (27) and (38). For this restricted 
domain the solution of (i) is therefore unique. 


2 • 2. Solution of the non-linear hyperbolic 
equation. 


In this section we consider the non-linear 
hyperbolic equation 

in the domain 

o < x < jt, o < / ^ T. (2) 


p(x) is > I, and p(x) and dpjdx are continuous and 
uniformly bounded in the whole interval, r is any 

given positive integer, while p^x, t), ^ 

are continuous and uniformly bounded in the 
domain (2). 


We wish to determine a solution of (i) which 
is regular in the domain (2), and which satisfies the 
boundary conditions : 

u(Oy f) ~ u(p, /) = o for all / in (2) ; 


(3) 
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The functions /i(>r), can be* expanded in 
a series of Sturm-Liouville eigenfunctions <P„ (X): 


Ux) - 2 Uix) = 2 (5) 

n n 

where we assume that the series 

2 X„li-„| and 2 (6) 

n n 

are convergent. 

For the solution of (i) we write as usual. 

u{x, /) = 2 »'„(/) <P„(x) (7) 

n 

The conditions (4) will be satisfied if the coefficients 
(/) satisfy the relations 

= ( 8 ) 

(« =- I, 2, ... CO). 

Further we have 

A (•''<■. 0 2 (0 ( 9 ) 

n 

where 

J?:« (0 "= 1 Pr 0 «’■ (^. ‘Pn dx 


2 { t , A . • • • > «'*1 (0 »•'*. (0 • • • ( 0 . 

...» 

( 10 ) 
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with 

(A {x, t) (x) ... 

(“) 

We assume, what we shall prove later, that the 
series 

^2 22 ^ 

2 K W <P« M and 2 , 2 ” q>„ (a') ... ( 12 ) 

n n 

arc absolutely and uniformly convergent in the 
domain (2). Then substituting (7) and (9) in (i) 
we obtain 

-^K (0 (>^) - 2 9 „ (.V) 

n n 

= Z :?:„(/) 9„ (.v), (13) 

n 

so that for all « > i we have 

- Q) == - (0 (14) 


The solution of this equation which satisfies 
the initial conditions (8) is given by 

(f) ^ cos V ^ sin V^„ t 

- sin V K - j) O) 4y 
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= <-„cos V^«^+ 

X 2 n> (j) ... {j) dy, (15) 

ky ... kr 

{n — • I, 2 , ... 00). 

This is an infinite system of non-linear integral 
equations for determining the coefficients a'„ (/). 
We transform the system by writing 

(^) = <■« cos V / -1- y” sin V ^ W. 

....”.. (16) 

/n if) = K COS V A„ / -f V K dn V K (l?) 

Then we get 

»„ (0 == - V ( sin V it -j) 2 

{/*i (j) + (j)} - I A, O') + O)} ••• (18) 

(« = I, 2, ... 00). 

We write for all « > i : 

gn Oj) = - sin V K it -y\ ( 19 ) 
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Then (ft) becomes 

= r gn ^ K ( Aj (j) + % O')} 

•'O 

•••{A,0) + ^*vO)}^. 

in = I, 2, ... oo). 

From (6) and (17) we get for all t in (2), 

S \ fn (0 I < ; (s absolute constant) (22) 

n 

Also for all « > I and all t in (2) we have 

I gn Oj) I = I - sin it -j) I < I (23) 

Further, from (8) § 1 . i we obtain (since x = o) : 
'Z\h„\< A h^pir), for all t in (2), (24) 

n 

and all k-^y ..., > i. 

From §2.1 we see therefore that one and only 


one solution of the system (21) exists for a restricted 
interval of /, and also that 2 | »„ if) ( is uniformly 

n 

convergent in this interval. 

Having found we get the coefficients 

wjf) from the relation (16). When these values 
of n>Jf) are substituted in as given by (10), 

we find that the series 2 | ^„(/) 1 is uniformly 

n 

convergent for / in o < / < T. Further, we 
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find from (i6), on account of (17) ind (az), tb^ 
2 "Kn [ uniformly convergent. 

n 

functions 'P„(x) are uniformly bounded for au 
xin o < < re and for all « ^ i. Therefore 

2 K ^n(.^) ^n(^) 2 :^„(/) (pjx) are absolutely 

n n 

and uniformly convergent in the domain (2) ; 
consequently from (13) we see that 2 -r.^- is 

also absolutely and uniformly convergent in (2). 


Thus 

u{x, /) = 2 »'„(/) q>„(.v) (25) 

n 


is the unique solution of (i) which satisfies the 
boundary conditions (3) and (4). 


Obviously there is no difficulty in extending 
the result to the equation 

dxK'' 'd.v) d/2 

where j is any finite integer. 


2 t) w'-, (26) 


? '^l 


2*3. Non-vanishing boundary values. 

In the theory developed so far we have as- 
sumed the boundary values to be zero. This was 
not done from any necessity, but only with a view 
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implicity* We shall prove now that non- 
shing boundary values are also amenable to 
treatment. 

Let it be required to determine the regular 
solution of 

5;KMax)-aT^=A(^>')“' W 

in the domain 

o < X < Jt, o / 5^ T, (2) 

satisfying the boundary conditions 

«(o, t) -= /) = •>1>2(0 all / in (2), (3) 

and the initial conditions 

u(x, o) = du(x, 6 )/df — f2{x) for all .v in (2). (4) 

Let fi{x) and f.Jix) be capable of expansion 
in Fourier scries 


/i(x) - 4 cpXv), (5) 

n n 

such that the series 

[< I and 2 Vx„ [<' ( 

n S 

are convergent. (6) 

In order that the conditions (3) and (4) are 
compatible, we must have 


H'i(o) 








O, 


( 7 ) 
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dt dfi °' 

Let^(x) be any function defined in o < >r < - 

1-1. d-g dh d^g d'^g 

such that & M= continuous. 


i(o) = |(o)=^|(o)=j^(o)=^J(o) = o' 

cf-\=zl \ = 

3 / ^ dx^ ^ ^ dx"^ ^ 3 ^ dx^ ^ 3 ^ 

d"^^ / 7t\ 

Then we define a function F(x, f) as follows ; 




F(x.ty 


P 2 {^> ^)=^2W+ 


for o < X < 
(x—Kf(d^% .,d% 




for X < Jt, 
3 


FA^> ^)+ {Fz(^, f)-F\{x, /)} (x- p 


r n zn 

tor < X < “ 
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Then we see that F, < - , -v-k and their first four 
d/ d/2 

derivatives w. r. to at are continuous. Also 

F(p, t) = \l)i(/), F(n, t) = %{t) (lo) 

dF 

The functions F{x, o) and . - (x, o) can be expanded 

U/ 


dF 


in Fourier series 
F(x, o) = S i; 9.(X). - (a-, o) = S F>.M, (ii) 

n ^ ^ n 

with 

2 ( i’' j and ^ [ convergent. (12) 


We define another function H(x, /) by 

, d2F 

I 


') = - a3 + +A(-- ') (■’) 

Then H, and t ^ are continuous in the 

dAT dx^ 

domain (2). It can also be verified that 

H(o, /) = H(a, /) = o (14) 


Hence H(x, /) can be expanded in a series 

H(x, /) = :^ /&„(/) <p„(a-), (15) 

n 

where 

2 V^„ I /^„(/) I is uniformly convergent. (16) 



62 


HYPERBOLIC EQUATIONS I 


Now we consider the function 

ip{x, t) == u{x, t) — F{x, t) ( 17 ) 

It satisfies the differential equation 



■-= A(-=^> ^)} + t) (18) 

We substitute for «(>:, t) in the right hand side 
of (18) by writing u — iv F, and get 

d ( , , , dw) ^hp 

— H (x, /) + t) T-V -f q'f (x, f)yp‘^ - \- ... 

+ (-v, t) 

= H (x, t) + 2 (.V, /) (x, /), (19) 

fc=l 

where the functions (x, t) can be easily calculated 
in terms of r, (x, t) and the powers of F (x, /). 
For example 

^(1) (x, /) = rp^ {x, t) F^-^ (x, /), 

?“ (^. 0 = *-7-' V. (■>'. 0 r'-“ (^. /). - (20) 


etc. 
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The funAion w (x, /) satisfies the boundary 
conditions 

ly (p, /) = » (o, /) — F (o, f) 

= % (/) - ■‘I’l (f) -= O, 

> ( 21 ) 

jy (jt, t') = u (jt, /) — F (a, /) 

= ■'I’j (^) — ^2 0 = ° 5 , 


u> (x, o) = « (x, o) — F (x, o) 

= /i M o) = W (say), 

o) = a" (^. o) - (-v, o) 




aF 

a/ 


a/ 

(x, o) ~ /j.^ (x) (say). 


i ( 22 ) 


We see from (5) and (ii) that the boundary 
functions (x) and (^) can be expanded in series 

(x) = Z //„ <P„ (x), h, (x) == hZ <P„(x), ... (23) 

n n 

and that on account of (6) and (12) the series 
^ 1 1 and ^ V I h" are convergent. ... (24) 

n n 

As shown in (26) § 2.2, the equation (19) for 
n>{x, t) can be solved for the boundary values (21) 
and (22), and then the solution u (x, t) of the equa- 
tion (i), satisfying the boundary conditions (3) 
and (4) can be obtained by taking 

u(x, t) = jp (x, /) -f- F (x, /) (25) 
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NON-LINEAR ELLIPTIC EQUATIONS 

3*1. Resume of the results for linear 
elliptic equations 


3 ‘i (i). The boundary problem for = o. 
Let D be a region bounded by a simple closed 
contour C, and consider the equation 






ay-°- 


(0 


It is proved in the potential theory that we 
can find a solution u(x, j) of (i) which is regular 
in D and continuous on C, and which is such 


that the expression 

au p 






takes given values on C. The curve C, the func- 
tion a, and the boundary values can have a 
fairly general character.* 


*r./. The articles on potential theory by Lichtenstein 
in the Enzyklopaedie der Mathematirchen Wissenscheften 
Bd. II (Analysis) C 
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yi (2). The Green’s function. 

The Green’s function for C is defined as a 
harmonic function (solution of \/^u = 0) of ■q 
which vanishes on C, and which is continuous 
in D except at the point (x,j) where it becomes 
infinite as — ^ log {(| — x)® + (t] — j)*}. 

Let g(x,j) be the harmonic function regular 
in D, continuous on C, and taking on C the same 
values as \ log j(5 — + (t^ — jY) . From the 
previous paragraph we know that such a function 
exists. The Green’s function is therefore 


G(l, Ti ; x,y) = g{x,y) — Uog{(5 — 

-f- (n -yf\ • (i) 

The function y) has derivatives t ^ and 

perfectly determined in D, and ^ are also 
determined in D except at the point 5 = 1I =y 
where they become infinite as ^zs^rTTr S2 


n —J 

(5 — xf + (11 —yf 


- (I _ ^)2 + _^)2 

respectively. 


Let 

^(•^> j) = I ^ 
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Then A(x,j) is shown to have ?. maximum M 
when (x,y) moves in D, and that 

o when D -> o. (3) 

Similarly let 

i 

Then it is proved that B(x, j) and C(x, j) 


have maxima such that 

I B(x, j) I < N, I C(x, j) I <5 N, (;) 

and also that 

N o when D ^ o. (6) 


Thus we see that M and N depend only on the 
curve C, varying continuously with C and tending 
to zero with C. 

3-1 (3). The boundary problem for 

= f(x, y). 

Consider the equation 

l\^u=f{x,j), (. 1 ; 

and suppose that we wish to determine a solution 
which is continuous in D and which vanishes on 
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C. It is easiiy seen that the solution is given by 
the integral 

u(x, j) = — ^ J /(I, n) G( 5 , d\dr\. (2) 

Suppose F is the maximum of /(x, j') in D 
and on C. Then we get for all x, y to. D and on C : 

lu(x,y)j<FM. (3) 


Further we have 

Then from (5) § 3 • i (2) we obtain for all in 
D and on C 

la"l<FN,| 5 ''|<FN. 0 ) 

The derivatives and ^ of (2) exist provided 

f(x, y) is continuous. But the continuity of /(x, y') 
is not sufficient for the establishment of the existence 

of the second derivatives - s and v“v> , and of the 

6x2 ^ 2 » 

relation /\^u — f. This will be the case if / satisfies 
the relation 

1/(^1. Ji) - /(^o» Jo) I < -K (6) 
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where (xi, and (atq, Jq) are any two points in D, 
r the distance between them, and K and two 
constants. 

However, if and are continuous in D, 
dx djy 

then u(x,y) has second derivatives which are them- 
selves continuous in D, and which remain continu- 
ous if the point {x,y) approaches indefinitely near 
to the contour C. Further, if 

... (7) 

then it is found that 



( 8 ) 


where A, and M, are two positive quantities indepen- 
dent of /(x, j), and where 

M, o when D o. (9) 


3-2. 


The non-linear equation 

T»/ 

= P(x>y,u,s^. 


dy/’ 


In this section we wish to determine a solution 
u{x, j) of the equation 
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2/7 ==: 


d^J( d^u 


P(x, y,u,p,q) (i) 



which is regular in the domain D 


and takes given values on the curve C. Sup- 
pose that the function P remains finite and well 
determined when 


l»|. 1^1. w 

for all x,j in D, L being a positive constant. Then 
we see that P will have a certain maximum jt. 

We now take the following successive approxi- 
mations in order to solve the problem. 

Take an}’^ function »i(x, y) such that 

Kl. IPi I. I < -L, (3) 

^ du„ , . 

^ dj 


and write 

J. Pv^i) (4) 

== P(>^, J, »2> p2. ^2) (5) 

AX = (6) 


Since we know the function »i(x, y), the right 
hand member of A) is a known function of x, j», 
say Thus (4) is of the form »2 = f(x,j) 
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which has been considered in § 3 * We can 

therefore find a solution u^{x, j) which is regular 
in D and takes given values on C. Then we can 
determine from the equation (5), and in 

general we can find Un {x,y) satisfying the equation 
(6) and taking given values on C. 

To establish the convergence of these ap- 
proximations we proceed as follows. 

Taking M and N to be the quantities defined 
in § 3 • I (2), we know that M and N tend to zero 
continuously with D. We choose D sufficiently 
small so that 

\iM<L,iiN<L. (7) 

Now P(x,y, /q, pj, ^,) has the maximum ji, so 
that on account of (5) and (3) § 3 ‘1(3) we get 

[ 1 < jiM < L, I I < liN < L, I ^2 I < H-ZV < I- 

But from (2) we know that as soon as P 2 
and ^2 between — L and + 1., the function 
P(x,j, p2, qP) will have the maximum ji. Then 
from the equation (5) we conclude that 
I «3 ) < < L, |y>3 I < |.N < L, [ ^3 ( < (xAT < L. 

Repeating this argument we see that in general 

\u„ \ <\kM < L, I < [xN < L, 

k„ ( < :'-N < -L. 


( 8 ) 
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(« I, 2, 3, 

We write 

3,4, •••). (9) 

and get from (6) and (9) 

— (10) 

V 2 = P(x,j, q„-^) 

- P(x,j, A- 2 ^n-2) (”) 


From (9) we see that all (x,j) vanish on C. 

We assume that the function P(Xyj, u, p, q) 
satisfies Lipschitz Condition in u, p, q, so that 

( P{^,y, th^Pi. qd - PQ<yy> "i>A. ^i) I 

<^|//2-«il + B|A-A + <^k2-^i|. ( 12 ) 

where A, D, C are positive constants. 


Thus, for all « we obtain from (ii) 
\V^Pn\<yi\^n~l\ 

+ + 

(« = 4 , 3, ...). 

On account of (8) and (12) we see that the right 
hand member of (10) has a maximum F such that 

F<iL{A + B+C) ( 14 ) 
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Hence from (5) and (3) § 3 ‘ i (3) we<»obtain 

<NF, llfKNF. 

Consequently 

-4 ( -s I + B I £ I + c I 1 < 

<RF, 


where 

R AM -h BN + CN. (13) 


From this we get 

I 1 < MRF, 1 1 <N.RF, I I < N. RF, 
and therefore 

-4 I ^ 4 1 + B I 1 + I I < 

Thus in general we find for all > 4 : 


rg\<N.Rr>-^F (16) 

From (16) we conclude that v„, will tend 

to zero, and consequently will tend to a limit if 
R = AM+BN+CN<i (17) 
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Now M 4nd N tend to zero with D, so that 
the Condition (17) will be certainly satisfied if D 
is taken sufficiently small. In this case the limit 
function u(x, y) is represented by the infinite series 

* — ^2 “b *^3 4 " ^4 + ••• + "b 


This will be a continuous function of .v, 
which has the first derivatives given by the series 


da _ da, , d2;3 

dx ~ dx dx + ••• + dx + 

, a^3 4- , 

dy ^ ^ ^ ^ 


1 

j 


(19) 


The function u (x, j) will take on C the same 
values as and therefore the given values. 


Further, since 

G(S, ii ; -v, y) d% di\y 


where G(§, t]; x,y) is the Green’s function for C, 
we get on passing to the limit n 00 : 

u(^,.y) - - P{1^ n, »,ll , 1 ^) 

G(i, i];x, y) dl d^ (20) 
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We have seen that a(x, j) has fivst derivatives 

^ ^ . We must prove now that the second 

derivatives also exist, and satisfy the differential 
equation (i). 


We assume that P(x, j, u,p, q) has first deriva- 
tives with respect to all its five variables. If the 
initial function u-^x, j) is chosen such that it has 
continuous second derivatives, then we see that 

the function/(x, j) = P(x, jy, will 

have first derivatives continuous in D and on C, 
and therefore U2(x, j) will have second derivatives 
continuous in D, so that P{x, y, qP) will have 
first derivatives continuous in D. 


We then calculate from the equation 

«2. P-2, ^2) - »i, Pv ^i)- (21) 

We have just proved that the right hand 
member of (21) has first derivatives continuous in 
D. On account of § 3 ‘ 1(5) this is sufficient in 
order that may have second derivatives. 

Suppose F is the maximum modulus of the 
right hand member of (21), and Fj that of its first 
derivative, then from (3), (5) and (8) §3*1(3) we 
obtain 
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[“'s I < iWF. 1^1 < NF, I I < 1 F + Af.F,. 

(^ 2 ) 

We go over to which satisfies the equation 

V®t'4 = «3. Ai. ^ 3 ) - J. «2. A. ^ 2 ). (^ 3 ) 

where we know from (13) and (15) that the right 
hand member has the maximum modulus RF. 


Let 

= Ma=. 1 *^ 1 . ir Max | C = Max 

(m) 


Then 


{P(-"^''>J^> ^3»A3» ^ 3 ) P(-'"'-%J'> P'i> ^%)\ 1 

< A'NF + B\lF d- A/jFi) + C'(XF -f MpF^) 

< (X A" + NA') F + MiA"F\, (25) 


where 


A" = B' -H a. 


(26) 
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Thus from (23) we obtain t 

I r J < MKF, I 1 < NRF, I 2^* I < IRF 

+ M,(aF+bF^, 

where 

a=^lA" + NA\ b = M^A" (27) 

Proceeding thus step by step, we see that the 

maximum moduli of ^ will be obtained 

" ox ox^ 

from those of ^ successively by the linear 

substitution of 

RF for F, aF + bF^ for F^ (28) 

We must make the substitution n times if we 
wish to obtain the maximum moduli of (x,ji). 

Thus the maximum modulus of 

ox ox^ 

will be 
dx^ 

XR" F + Afi{<z(R“ + R"-^ b + R ”“2 b'^ 

+ F + b”"^^ Fi) (29) 

Now R and b contain M, N or Mj as factors, and 
each of these tends to zero with D. Thus taking 
D sufficiently small we can make both R and b 
less than unity, so that we find 
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dx^ 


-> o when ft 


00 . 


(30) 


This shows that 


d^u 

dx® 


exists and is continuous in 


Z 7 . Similar results hold for ^-5 . 

dj2 

We conclude therefore that the function 
given by (20) satisfies the given differential equation 
(i), and since it has been shown to satisfy the 
boundary conditions also, it is the required solu- 
tion of our problem, provided that the domain D 
is sufficiently small. 


3.3. Uniqueness of the solution. 

Suppose j) and j) are two solutions 


of the differential equation 

= P(x, y, u,p,q), (i) 

which are regular in D and take the same boundary 
values on C. Let 

v{x, j) u.lx, j) — Ui(x, y) (2) 


Then v{Xy y) satisfies the differential equation 
SJhf P{x,y, «2, A. ^2) — J. «i. A. a) 
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where u lies between j) and y). 

We see also that v is regular in D and vanishes 
identically on C. 

We shall prove now that v will vanish identi- 
cally in D provided D is sufficiently small, and that 
this zero solution is the only solution of the equa- 
tion (5) taking zero values on C. 


We assume that the function P and aU its 
derivatives occurring below are continuous. We 
define three functions d, e,f by 


j)= - h Qp , e{x, jy)- -k /(-V, j)=- - 

The equation (3) then becomes 


dP 


( 4 ) 




zd 


dp 

dx 


dp 

ze ^ ~r p 


dj 


o. 


(5) 


If p(x,j) is a solution of (5) in D, then evidently 


n tdV dp 


dp 


( 6 ) 

Integrating (6) by parts and remembering that 
V = o on C, we get 




dd de 

_i_ 


(7) 
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af \ I r 

_dP J a 2 p J sap 

S* ^ Qxbp ^ djdq’ 

If 6(>r, j/) is positive in D, i.e., if 
SP I S2p j sap 

S« ^ SA-Sp ^ ^ ° 


( 8 ) 

( 9 ) 


n D, then from (7) we deduce that for all x,j in D 


V = o. 




(10) 


This means that v{x,j^ vanishes identically in D, 
so that «2 j) = «, {x, f) in D. The uniqueness 
of the solution is therefore established for any 
arbitrary D. 

In the general case when the inequality (9) is 
not satisfied, we employ the following procedure 
to establish the uniqueness for sufficiently small D. 

Let j) and y) be two arbitrary real 
functions which are continuous in D and on C. 


Then since v — o on C, we have 

1 L (‘P + aj ^ °- 


(XI) 
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Adding (7) and (ii) we obtain 


IL{( 


2 Qp dj, 

> +(5,) + 


6 x 


~ dx ' ^ 


f^^^dxdy = o. (iz) 


The integrand is a quadratic form in v, ^ , 
and will be definite if 




(15) 


We shall prove that if 0 (x,y) is continuous, 
we can determine two functions 'P(x,jv) and j) 
such that in the neighbourhood of an arbitrary- 
point the relation 


(p2 -j- ip2 — 0 


c 5 <p 

dx 


dij) 

dj 


is satisfied. 


(14) 


To prove this suppose that r?P- is the maximum 
modulus of — 0. We wish to determine two 
functions <p and ^|) which satisfy 

(15) 

The inequality (15) implies necessarily the in- 
equality (14). We can easily find one pair q>, \l> 
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which satisfies»(i5). Thus put = o, and take q> 
a function of x only, satisfying the equation 



where m\ > The solution of (i6) is 

<p = tan {myx + a), (17) 

where a is a constant. 


Thus we can certainly determine 9 in an 
interval between two parallels to the j/-axis whose 

mutual distance is less than — , and this q> will 

my 

satisfy the required relation. For this pair q>, the 
inequality (15) and consequently the inequality 
(13) will be satisfied. 


Hence, in a region of the plane round the point 
for which we can determine the continuous 
functions <p, ip satisfying the inequality (15), we see 
that the equation (12) will be possible if and only if 





(18) 


This shows that provided D is sufficiendy 
small, the differential equation (3) can have no 
solution, other than zero, which is regular in D 
and ''iiJiich vanishes on C. 
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Consequently, there is one ?nd only, one 
solution of (i) which is regular in D and which 
takes given boundary values on C, provided D is 
sufficiently small. 

The restriction that D should be small, can 
be removed if, instead of the general equation 
(i), we consider the particular non-linear equation 


= P(^y J, «), 


ap 

du 


> o. 


(19) 


* It can be shown* that for any arbitrarily 
large domain D, the equation (19) has one and only 
one solution which is regular in D, and which 
takes given values on C. 


*P. Frank u. R. von Mises: Die Differential-und Integral 
gbichungen der Mechanik und Physik. vol. I, Chapter i8, 
§ 2 and § 3. Also E. Picard : Journal de Math. (4 series) 
Vol. 6 (1890), Chap. Ill, p. 173. 



CHAPTER IV 


NON-LINEAR PARABOLIC EQUATIONS II 


4*1. Resume of the results for linear 
parabolic equations. 


4’i (i). Fundamental solution. 

Consider the equation 
62 « du 


(i) 


in the domain D bounded by a portion A^ of 
the x-axis and the two curves and Cg given by 
the equations 

Cj : X h-^ (/), (^2 • ‘V = ^2 (0> (^) 

and bounded above by any parallel to the 

x-axis. Without restricting the generality, we can 
assume that the time A^ A^ is given by / = o. 


The equation which is adjoint to (i) is given by 




( 3 ) 


Let P and ^ be any two points in D whose 


7 
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coordinates are (x, t) and (i, t)) resp/?ctively ; then 
the fundamental solution is 

U(l,y]; X, t) = -7-^- r (4) 

y t— r\ 

It can be easily verified that { 7 (?, ti; .v, /) is a 
solution of 5« = o considered as a function of x, /, 
and of Sji' = o considered as a function of 5, 4. 

Let x-^ and x^ be the abscissae of A-y and A 2 
respectively, and consider the integral 

fa;o T _ 

Kix, y a^, (5) 

where is a function continuous in (xj, .Vg). 

This integral is called ‘Poisson’s formula,’ 
and it represents a solution of the equation 6« = o, 
which is regular in D except on the characteristic 
Ay A^. It vanishes at all points below this charac- 
teristic, and also at aU points of this characteristic 
outside the segment Ay A2. In the portion of the 
plane situated above Ay A2 it is a holomorphic 
function of x, t. 

We have further 

lim K{x, /) = 2 y jt if •’^'o is in (-^i, 1 

(x,/)-»-(xo,o) i- 

= o, if Xq is outside (x^, x^. j 

(6) 
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4'i (2). boundary problem for 5u = o. 

With the help of a system of two linear integral 
equations it can be shown that one and only one 
solution of 8« = o exists which vanishes on A^, 
ind satisfies on Q, relations of the form 

h (0 57^ + Ifi (0 « (/), (i = I, 2) (i) 

It is necessary for the truth of this existence 
theorem that the curves Q and Cg are rectifiable, 
ind that their equations satisfy the inequalities 

1-f-a 

I 1 < -t- 1 ^ — n h , (o < « < i) ( 2 ) 

where L and a are two constants, and / = i, 2. 

The following two properties of this solution 
of = o are analogous to those for the potential 
functions. 

1. If a solution is regular in D and continuous 
on (C), it cannot surpass its extreme values which 
it attains on (C). 

2. If a sequence of solutions converges uni- 
formly on (C), it will converge uniformly in D 
and will represent a solution. 

4*1 (3). Green’s function. 

Through the point P with coordinates (x, t) 
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draw a line parallel to the Ar-axis cutting the curves 
Cj and C2 at and Mg respectively. Let the 
domain bounded by M^ A-^ Mg M^ be called D^. 

Suppose H(i, T] ; x, t) is a solution of the 
equation 




(0 


which is regular in D^, vanishes at each point of the 
segment Mj Mg and takes the same values as the 

I __ (r-C)2 

function rjj- — e on the curves Af, A-, and 

\ t — T] ^ ^ 

•Afg^/^g. 

Then 

1] ; I) ; ^) — n ; -v, (2) 

is called the “Green’s function” of 8« o for (C). 
It is the solution of hu = o in a-, t, and of 81 — o 
in I, t]. It vanishes when the point (i, n) [or the 
point is on Q or Cg, and behaves in the neigh- 

bourhood of the point (a*, /) [or of (i, tj)] as the 
fundamental solution U. It is always positive in D. 

It is easily seen that the Green’s function is 
symmetrical in (a-, t) and (5 , ti) : 

G(x, / ; I, Ti) = G( 5 , n (3) 

Further it can be shown that 

G=g Ti ; X, /), (o < ^ < i) 


(4) 
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where ^ is a Jsounded function lying between zero 
and I. 

Again 

OG _ 5L7 dH , . 

dx Qx dx ’ 

from which we see that = o when (5, t]) is on 

Gj or Gg and (.v, /) is in D. Similarly all derivatives 
of G with respect to a: and / vanish when (§, t^) is 
on Cj or Cg and (.v, /) is in D. Conversely, all 
derivatives of G with respect to ^ and ri vanish 
when (x, /) is on or and (|, ri) is in D. 


4 ’i ( 4 ). The fundamental formula. 

Let (p(.v, /) and 4’(-v, /) be t w o functions having 
second derivatives. Then we have 


wl»69 — <P8pl' -- 1' ( 


0‘'qp 


Oa* 


( Ip 


oqp 

6/ 


) - 


82\p 

5x2 



5q) 

5.V 



5 

5/ 


(<pip). 


Hence from Gauss’ theorem we get 



— 98 in') = 1 99 t/.v + ( '9 

J «;) ^ 


O 9 

dx 



(I) 
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This is called the extended Greek’s formula, or 
Riemann’s formula. 

From (i) we get, if u(x, /) is a solution of the 
equation 8 » — f(x, /), 

(a) 2 V A u(x, f) 1 

(P) 

(y) o j 

\u{\, ii) d\ 

“ J I/) V/^__ *l) W 

This is the “fundamental formula” in which (a) 
holds when {x, f) is inside D, (P) holds when (,v, /) 
is on (C), and (y) holds when (.v, /) is outside D*, 

4‘i (5). The first boundary problem for 
8 u = o. 

We wish to determine a solution of 8 « = o 
which is regular in D, and which takes the following 
boundary values on (C) : 

a ^ tPj(/) on Cl, u ^ <P2(/) on C^, (i) 

* For details see Goursat, Course d’ Analyse, Vol. 3 . 
§ 544 PP- 508 et seq. 
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and 


u — (p(x) on (a) 



where and are two numbers such that 
< Ag < tvhere tp(a-) is a function continuous 

in (<3i, ^2) coinciding with qj(.v) in (a-j, a- 2). From 
§ 4’ I (i) we see that 6« — o in D, and that u — <p(x) 
on ^1^2- 

Suppose «(a-, /) is a solution of the problem 


(i), (2), and write 

»(x, t) ^ 'u{x, t) -I- /) (4) 

Then Wj (a-, t) satisfies the equation 

8«i ^ o, ( 5 ) 
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and the boundary conditions 

% == <Pi W on Q, = q>2 (/) on C2 ; ( 6 ) 

Ui = o on ( 7 ) 


provided of course that <Pi (o) = q>2 (o) = o. 

To find we have from the fundamental 
formula (za) § 4 • 1 (4), on writing f — o, 

( t -£)2 
4 («— 77) 

{^1 + {%- - " afj- V } 

Similarly, if in (i) § 4*1 (4) we substitute 
( 5 > 'h) fo^ 'P iind H (^,r\; x, /) for ip, and remark 
that «i = o on and H - - o on Mj M^, we get 

o = — [ U-, H 

+ ■’")</, (9) 

Subtracting (9) from (8), and remembering 
that Hhas the same values as U on Q and C2, we get 


»i t) 
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The values of are given on the curves 
Q, Cj, so that the integral (10) gives us (x, /) in 
D. Hence on account of (4) the required solution 
u (x, t) is obtained. 

Any function /(x) is said to admit an incre- 
ment of order y (non-zero) if for any two values 
Xj and Xg (> Xj) of x, we have 

1/(^2) — I < L I X3 — Xi (ii) 

where L and y are constants, and y is positive. 
f(x) is called by .some writers to be continuous (L). 

Thus if qp;(/), <p'(x), ^ 4(0 

exist and admit increments of order non-zero at 

all points of D and (C), then and exist, 

and we can determine a number y < i such that 
for a given increment / of / we obtain 

|A«!<L|a/|. lAs^l<L|A<r"", 

!Aj"|<L|A«|’' (12) 

where L is a constant, and where is the smaller 
of the two numbers I and i — y. 

All the results of this section hold also for the 
equation 8« = /(x) 
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4.1 (6). The equation 6u = fi(x, t). 

We consider now the equation 
^ dhi du j.. 


(0 


and determine the solution u{x, t) w'hich is regular 
in D, and takes the same boundary values as (i) 
(2) of § 4.1 (5). Then by a method similar to that 
employed in the previous section we get for the 
solution of this boundary problem: 


u{x, t) — 


7 ~ 7 ~(, "(§> n) 


o(^(s» 7j; Xf t) 


dy\ 


+ IT/ a 

2 V TtJA^A, 

— ^ I n; -V, /)/(?, T))^/i di]. (2) 

If the function f(x, /) has the maximum modu- 
lus F in D, and if tpC-v), ‘P'(.v) have maximum moduli 
M, M' respectively in (xq, a 2), then we find in D 

\«\<M+Ff,\ll\<M'+ ( 3 ) 

Butif l/i < F (F constant) then we obtain 

|«|< T" I < Af' -h L' /». 

(o < j8 < i), 


( 4 ) 
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where U is a plumber depending on the increments 
of (/), q>' (x) and b. (/), and is smaller of the 
two numbers ^ and i — y. 

Now suppose that (x, /) is a solution of 
h u —f vanishing on (C), so that on account of 
(2) we have 

«0 (a% ^ ^)/(^’ 

(5) 

Then it is found that if \f\<,FPm D. 

I "0 1 < F . I a"" 1 < Lra (i, V + I) 

( 6 ) 

where B denotes the Beta function. 

If ( f(x, t) satisfies the condition 

i A/i --■ l/(-v, t -i- A t) —/(-V, 

f{x, o) -= o, (7) 

K being a constant, then we obtain 

lA«ol <-LKAf,:Aa 5 'l <LK(Aty^‘,' 

I A ! < LK (A /)’. I (8) 

|«„| <LKr*\ 1 <LKr, 

I u 1 >1 I ’ ' Oy ' 


where A is a constant. 
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4*2, The non-linear equation 8 ,u=P(x, t, u, p) 
4*2 (i). Existence theorem. 

Consider the equation 


bu 




e/ 


P(x,t, u,p\[p^l~] (l) 


where the function P remains finite and well deter- 
mined when I A I < iV and \ p\ < N, provided 
{x, t) remains in the domain D. P is a continuous 
function of a-, /, satisfying the conditions 


I PCI, ri)-P(|,/)l<i:(/- 
|P(|,q)-P(A- ,^)\<K\x- 


nl'M 


i'V 


J 


^(o < Y<l) (2) 


P satisfies also the Lipschitz condition in «, p, 

i.e., 

I -PC-v, /, A) — 

< A I Ag i ■ r ^ I p 2 A I (3) 

We wish to determine a solution u (x, /) of 
(i) which is regular in D and which satisfies the 
following boundary conditions: 

a = (Pi (/) on Cj and » = tpg (/) on Cg, (4) 

and 

u (x, o) =: q) (x), i.e., a = tp (a-) on Ay A^. (5) 
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The curves Q and Cg are such that their equa- 
tions satisfy the inequalities (2) § 4’i (2). <p'(x) 

exists and admits of an increment of order non-zero; 
q>j (/) and qjg (0 ^^dmit of increments of order > 4, 

Let Aq (x, /) be any arbitrary function such 
that I Aq I < N and | />„ 1 < N, where we write 

p„ — (ti = o, I, 2, 3, . . . . Then determine 

successively the, solution of the following equations: 


5«i=--P(x, A o,A). (6) 

(7) 

8»„ = P(x, /, A-i). (8) 


Ag, ••• taking the given values (4), (5) on 

(C). 

Since we know Aq (a*, /), the right hand member 
of (6) is a known function of (x, t). From 4‘i(6) 
we can therefore find a solution of (6) satisfying 
the boundary conditions. Moreover, since | Aq |, 
I A 1 < the function P(x, t, Aq, A) has a maxi- 
mum modulus F. On account of (3) and (4) 
§ 4 • I (6) we get therefore 

1 a, 1 <M+F/, |Ai<^P-FL'/^ (9) 

Af, M\ \J jS having the same significance as in 

§ 4-1 ( 6 ). 
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Now, «i and must not lie outside the interval 
( — N, N), so that we must have 



/ must be therefore less than the smaller of these 
two numbers if | \ is to be less than N. for this 

sufficiently small domain wc then solve the equation 
(7) and get the same inequalities (9) for A’ which 
on account of (10) gives us again j % ] < N, \ p^l 
< N. Similarly for and other approximations. 
Thus we are assured that for the sufficiently small 
domain (10) //„, 7),, remain inside the interval 
(- N, N). 

Now we write 


v„ {x, t) = 

= »n 11 fo, /) - 


0, 

(« 

-- I, 2 , 5 , . 








(XX) 

and get from (8) and (i 

X): 





8/^1 = 

P{x, t, u^,pP) - 

— P(x. 

.A 


Po) 

(X2) 


P(X, /, «2, A) ' 

- P{^’, 

, A 


A) 

(X3) 

hv„ = 

P(x, /, u„,p„) 

- Pfo 

•.A 

A 


(X4) 


Evidendy all {x, f) vanish on C. 
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Since P sjitisfies the Lipschitz condition (3), 
we get 

I 1 < ^ i I + C 1 A - A I. (^5) 

\iH\<A\v,\ + C\^^\, (16) 


1 5'’, I < ^ ! ".-1 H- c I 1 (. 7 ) 

Since | i, |/> |o, i »i |, |A I < we find 
from (15) that the right hand member of (12) has 
the maximum modulus Mq such that 

M^<zN{A\-C) (18) 


Then we get from (6) § 4’i (6) on writing y — " o 
and F — Mq : 

I r, I, I i < L B (J. .) fll- 


< 


Zl^ ^ 


Ml (say) (19) 


From (16) we see that the right hand member 
of (13) has a maximum modulus Afj (A -j- C), 
and therefore for the solution we find: 


1 1. I S I < + C). Ji(J, I + i) 


r(2) 


(20) 
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Proceeding thus step by step, we get in general 

\ 2 ' 

We assume the domain to be so small that 
L (A + C) V nJ < i, i.e. / < 

dp 

Then v„, will tend uniformly to zero, and 

consequendy will converge uniformly to a limit 
u (x, /) represented by the infinite series; 

u (Xf /) = ~h % h ••• 4 ” + ••• (23) 

This function u (x, t) is then the required solu- 
don of our problem. 

4*2 (2) Uniqueness of the solution. 
Suppose the equation 

hu P{x, t, u, p) (i) 

has two solutions u■^ {x, t) and u.^ {x, t) which are 

regular in D and take the same boundary values 

(4), (5) of § 4.2(1) on (C). We write 

V (x, t) = ix, t) — (a-, /) (2) 

and see that p is regular in D, vanishes identically 
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on (C), and satisfies the differential equation 
8 t> = P( X, /, U 2 ,p 2 ) — P(^> A 

= (say). (3) 

Let M be the maximum modulus of g (x; /) 
in D, then we get from (6) § 4.1(6) : 

IHand! J^|<LMB(i,i)/J/2 ( 4 ) 


But since P is Lipschitzian, we get 


l^(x, /)| <^ + C!^^; ( 5 ) 

<(A-h C). L MB (4, i) p/2 

Ml (say). (6) 


Since the right hand member of (5) has now 
the maximum modulus M,, we get 

IM, <LAfjB(S,i+ 


(A + %) r(2 ) {l(. 4 + C) V n/j‘= A4(say) (7) 


V 1 and 


Proceeding thus step by step, we get in general 
Si' , ^ M 

^ (^ + 01(1 + 1 ^)' 

{L -j- C) jt/}**, (8) 


8 
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But L(A + Q V Jt/ < I on aacount of (22) 
§4.2 (i). Therefore we see that for all x, / in D 
the function p is identically zero, and consequently 
the two solutions and Wg ^re identical. 


4 ' 3 . The non-linear equation 8u — P(x, t, 

u, p, q). 

We consider the equation 

Su=P (x, /, «, p, q\ [ , q - J"] , (0 


for the boundary values 

(/) on Cj, u <P2 (/) on Cg (2) 

and 

u (jv, o) = qp (pc') for X in (3) 


We assume that the function P (x, t, p, q) is 
continuous for x, / in D and when | » |, |p |, | ^ | 
< N. Under these conditions the four derivatives 


dP 6P 8P dP 
hi ^ du Qp ’’ ^q 


We assume further that 


these four derivatives satisfy the Lipschitz condi- 


tion, and that 
ed, so that 


6P 

dq 


is essentially positive and bound- 


ap 

dq 


< K- 


(4) 
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The curvos Q and Cg ate such that h\ (/) and 
if) exist and admit increments of order non- 2 ero. 

Similarly, <Pi if), 92 if)-, 9' {x\ 9" (x) exist and 
admit increments of order non-zero. Moreover, 

I ‘Pi W > ! ‘P 2 W I. I 9 (at) !, 1 9' {x) I are always less 
than N. At the points (>rj, o) and (Xg, o), 
we must have 

9" (.Vj) - 9'i (o) - 9' (.Vj) h\ (o) 

= P {.Vj, o, 9 (-Vj), 9' (a'i), 9'i (o) -- 9' (aTj) h\ (o)}, 

‘P"(-V2)-9Uo)-9'(Xg)//(o) 

= P {Xg, o, 9 (-Vg), 9' (Xg), 92 (o) — 9' (Xg) //^ (o)} . 

(5) 

Finally, the equation 

^>" (x) — 0 = P {x, o, 9 (x), 9' (x), 0] ( 6 ) 

should define a function 0 (x) lying between — N 
and + N, and admitting an increment of order 
non-zero. 

Wc solve the equation (i) again by the method 
of successive approximations, and start with a 
function /q (xq t) satisfying the equation 

6 «i = P {x, o, 9 (.v), 9' (x), 0 (.v)} , ( 7 ) 

and taking the boundary values ( 2 ), ( 3 ), on (C). 
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For further approximations we sep> if 






dx 




5/ 


bVi -- P (x, /, qP) — P (x, /, (p, q)', 0), 

»2 = «i + (8) 

6t'2 = P (:xr, A »2.p2. ^ 2 ) - -P A "i,Pi, ^'i), 

*3 = + ^2> (9) 


62^„ = P (a-, t, />„, ^„) - P (x, /, «„_i, 

Pn-1» ^n-l)> ^n+1 ®n "i” ^n> (^°) 

all the vanishing identically on (C). 

We form the function 

«« = «i + d- H H- ••• T- . (ii) 
and see that »„ will be the solution of (i) taking the 
given values on (C). 

Now from ( 10 ) we obtain 

= -P ^»-l + Pn-H ■ 5^^ ?n-l + ^ 

- P(-^> P »«-l. P»-l. (12) 

so that enters on the right side. We must 

hv 

therefore establish the convergence of 2 
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for \\'hich we^ must find the maximum modulus 
of . We can do this by utilising the increment 
of the right hand member of (12). For this it is 
necessary to find the increment of • itself. 


We write for brevity 

a7^=-7„-i ( 13 ) 


We shall calculate the increment of rp„ when 
i undergoes an increment A A consequently 
when ... undergo the corresponding 

increments A A ••• • We shall have 
therefore to apply the formula for finite differences 
of a function of seven variables t,u,p,q,v, 
where we have suppressed the index (« — i). 
We denote by t', u', p\ q', v', w'y j' quantities 
which are contained between the extreme values 
of t, u, p, q,v,w,j, and write 

?<!' ^ -P{x,t',u',p\q'\ 

. - ip (a-, «' -1 v\ f q- / + /). 

We shall have then 


A = A ^ ( a/ 
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+ A/f' 


,dP‘2> 
'y~dp'^ 

dp<2» 

+ A - qy 


ap<i) 

dp 


ap<2) ap(i) 
ap( 2 ) 


(I) ,01''“' or'",. 


ap< 2 ) - 

f A»' -^+ Aj 'ay • ('4) 


We assume that the arguments of the function 
•»!> always remain in the field of variation which we 

ap 

have defined, so that under these conditions < - has 

dq 


the maximum modulus Mq, and 


6P dP 


have the 


maximum moduli M^. Since 


du ’ dp 

dP dP dP dP 
5/ ’ d// ’ 5p ’ dq 
are assumed to be Lipschitzian in »,p, q, therefore 
for an increment p,'u, /s'p, [\'q of «, p, q, the 
increments of the four derivatives of P are given by 


A I I A I 1 A I A 

^ d/ I’ I ^ du I’ I ^ dp dq 

<M2(|A'«1 -f- I A'PH- 1 AVI)- 


(15) 


Let \v], \w\, [ y] denote the maximum moduli 
of V, w,y in the interval A A then from ( 14 ) we get, 
on reintroducing the index « — i, 

I A 'Pn I < -^1 (1 A ^n-1 I + I A ^n-1 I + 1 AA-1 D 

+ M2 {[p,._i] + K-j] + {| A/ I 

+ I A ^«-i 1 + I A Pn-i I + I A 1} • ••• (i6) 
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We study now the behaviour of the succes- 
sive approximations, and start with the function 
»i (^, t) given by the equation (7), which is of the 
form hu — f(x). Thus from (12) § 4.1 (5) we get 
for > o : 

I A «1 1 < L A A I A s"'! < L (A 

|a|"'I<L(A0' (>7) 

Suppose the domain D is sufficiently small 
given by the values of / in o < / < /, then from 

(17) we get, since Uy (.v, o) = 0(x), 

! «i — «P(x) 1 < L/, I — <p'(x) ] < LA+^, 

1 ^ 1 - 0 (a-) I < LA. (18) 

Now we take the equation (8), viz. 5^^ = tJJi, 
and apply the formula (16) for « = i. We also 
write 

»o = ‘pM. a = % = ^ (•^). 

1^0 = «i - «^o A - A. 7 o = - ^0. 

so that 

A = A ffi, A»^o = A Pi, A Jo = A qi, (19) 

since «o>A>^o independent of t. From (16) 
for n — i and from (17) and (19) we get therefore 

lAAl<-?^i(A/r (20) 
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Consequently from (7) and (8) § 4. 1(6) we obtain 

1 A I < M' A I A 1 < M' -^1 (A 1 

I AjiKjtiCiCA/)". j 
I I < H iCi In'rK^iK, A+i/-'-, [ 

I Jl I < H- ^1 J 

where n depends on the given values and the Con- 
tour. We can always choose y so that it is the same 
in the formulae (17) and (21). 

Suppose that in general 4»„_j is zero on 
and that we have 

1 A n-n-i I < fC„_i (A /)". (22) 

and that analogous to (15) we have 

I A v„-i I < ^ K-i A A I A w'n-i I < ^ 

X (A 1 A j„-i I < K-i (A 
I I < ^ /-i. I I < A+v., i ^"5) 

1 Jn-l I ^n-1 A. j 

We can then write 

®n-l = "1 + + ••• + 2^„-2. 

^n-2 — ^] + f^2 + ••• -f- (^ 4 ) 

so that from (23) we get 

! »»-i — «i I < ^ i Pn-i — A ! < 

J* I I < (i i’„_2 A. (25) 



NON-LINEAR EQUATION 


107 


From (17) and (25) we find then 

I A ««-i I < (L + n A A 

1 A Pn-X I < (L + (A 

I A <ln-x 1 < (L + (A '^n-a) (A Pp ( 26 ) 

Consequentiy on account of (16) A '*J’« will 

be of the order y, and we can write 

\A%\<Kn(A tpy (27) 

where is obtained by the relation 

^ (A + -^fo -^4 A 

x”(i + V/+ /) l(A // + (A -f (2 A + i)} 

(28) 

Thus the inequalities (21) would hold for the 
index n with the same y, and if the expression (28) 
is less than i, the series ^ will be convergent, 

n 

and consequently the series Z v^, 2 will 

n n n 

also be convergent. 

Hence it is necessary that we should be able 
to determine an upper limit to the expression (28) 
which may be equal to q < i. 

In this case 


i’„_2<fCi(i + e4e*‘+... + e"-=*) 



(^9) 
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SO that on account of (i8), (24), (zj) we get 
I # — <P (x) I < I — <P (x) 1 + 1 + ... 1 


<L/ 4 -^ 1 -A+\ (30) 

and similarly 

\p-'t'(x)\<L /"» + P*<i% 

\f^e(x)\<LP+ A/y (,,) 


Substituting (29) in (28) and remarking that 
we shall get an upper limit for ^ , and our 

^n-X 

object would be realised if this limit is equal to q, 
for then we shall get successively „ " < q for all «. 

^n-l 

We obtain thus a relation of the form 
Q = + C /\{ q < x ), (52) 

where A, B and C are polynomials in /. 

Now (32) is an equation of the second degree, 
which gives us two roots between o and i if 
B/^ + a^ < I, and (B/o + C/^ - 1)2 

— 4^/^ >0. (53) 

This gives us a maximum value /j of /. 
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We must now write the condition that the 


values of shall lie between — N and -|- N. 

Let 

M = Max I qj (x) \, M’ = Max 1 9' (-v) |,, 

M" = Max d (x), (34) 


so that M, M' and M" < N. Then from (30) we 
have 

M+ L/-f --^-- 1 - < N, 

I 9 

M' -b L/^ f - ’‘^"1 ■ /" ^1/2 < 2v, 1 ... (J5) 

I — Q 

M" -f- L/y f < N. 

1 — Q • J 

We must therefore take a value of g such that 
it satisfies the equation (52) and the inequalities 
(35). Now these inequalities are satisfied by / = o, 
and Q = o, so that we shall certainly get in this 
way another upper limit /g. 

For the height / of our domain we then take 
the smaller of /j and 4. Then the solution would 
be determined in a sufficiently small domain of 
height /, because the series 
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will be absolutely and uniformly convergent only 
when / < /. The terms of the series are conti- 
nuous functions and therefore the sums will also 
be continuous. 

The solution of the problem is then given by 
u {x, t) = lim (at, t), (36) 

n -> 00 

where u^{Xyt) is given by the relation (ii). 

The uniqueness of the solution can also be 
established as before for a sufficiently small domain. 



CHAPTER V 

FURTHER RESULTS FOR HYPERBOLIC 
EQUATIONS 


5'i. An infinite system of non-linear In- 
tegro-differential Equations. 

Consider the system 

«« (•^) = (* in J) I'nijl »r (j) I < (A)} (0 

Jo 

(« =- I, z, ... oo), Avhere {y) denotes the derivative 

and ( j), «' {y) represent the totality of the 

infinite sequences. We stipulate that the solutions 
(x) as well as their derivatives /4 (•><■) should 
belong to Hilbert space, i.e., that 

n n 

should be convergent for all values of x in any 
finite interval. 


^in 

dx 


The functions g„ (x, y) and their derivatives 
are continuous for all x and jino<jy<x’<oo; 
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and for all values of x > o and all, real values of 
belonging to Hilbert space, the functions 
(j; a^; are uniformly continuous and satisfy 
Lipschitz condition; 

I f . ( j; “r, P.Oi <)) 

(2) 

/c=l Ar=l 

where and j 3 „ are positive constants. 

Differentiating both sides of (i) w. r. to x, 
we have 

< (•^) = gn ^'n K (•^) ; K ('■^0 1 

+ 1* ( 3 ) 

Two cases now arise accordingly as (x, x) 
does or does not vanish identically for every n. 

5 '2. The case when g„ (x, x) == o. 

Writing (x) --- /4 (x), i, 2, ... 00), we 
get from (i) and (5) §5.1 the doubly infinite 
system of integral equations 

»» W = 1 .?« y) Pn {yi »r (y ) ; k ( y)}4y^ (0 

K (^) = j* ^ y)F«{y>»r ( j) ; ( j)} (^) 

(« = I, 2, ... OO). 
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We solve ^this doubly infinite system by the 
usual method of successive approximations, and 
write for this purpose 

P„ { j; ^ 

Jo J 

and 

Pn\^) O. ] 

Pr”'~^Kj)}4y [ 

We shall prove that these approximations are 
uniformly convergent for all x, so that 

lim «<”•> (x) -- «„ (.v); lim (x) (x), ... (5) 

n -> 00 n -> 00 

and that 

2 ul (x) and 2 (a-) (6) 

n n 

are convergent; further 

lim 2 {«!’"’ (at) - //„(x)}2 

m -> 00 n 

-- lim 2 (.v) - (x)} 2 o, (7) 

w->- c» n 

uniformly tv. r. to x. 

To establish these results it is necessary to 
impose more stringent restrictions on the given 
functions than in Chapter TI. Thus we assume that 
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constants b, and v exist {a, ^ >.o; n, v > o,) 
such that 

2{[ I 

" / \o (8) 

and 

2 \<^^gn(g<,y) 4y\ 

2 \^^l\^gl{x,y) 4y\ 

Then from (3) we have 
«<,”+») (x) — (x) = j* (x, y) [F„ {j; «<”•> (j); 

— Pn {jl 

so that 

{ff(m+l) __ 2 

< ( 4y • r [“n 2 i«lr\y) - «*“■“ (jv)} ^ 

•'0 Jo k 

+ ^» 2 (j)} ^ (10) 

k 
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Summing over^ and taking account of (9) we obtain: 

2 {«(«+!) (x) — w'™' (x)}2 

n 

< bx " r 2 (7) - (7)} 2 

•^0 n 

+ (7) — (j)} (”) 

Similarly from (4) we get 
2 (x) — (.v)} 2 

n 

< [ 2 [ (7) — ( y)} 2 

Jo n 

+ (7) — (7)} (12) 

We write 

= 2 {ulr^^^{x) - /4”»(x)}2 

n 

+ 2 (13) 

n 

and obtain from (11) and (12) 

^m+l (x) < 2 /;x' r (7) (14) 

Jo 

Hence 

iPi (x) 2 (a-)} 2 + 2 (x)} 2 

n n 

^ 2 x**, ^ (8). 

Therefore from (14) we get 

7^2 (x) < zbx" [ zaj'^ dj — 


9 


4^^x“+’'+i 
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Similarly 


n^^(x) < zbx^ [ — dy 

•^0 M- I ^ 

^ Sab^ x*+2’'+2 _ za . (zby . 

■ 11 + ^+ 2 ~(|A+ i)(m.+ I + v + i) ■ 

In general, we get by induction for all > i ; 
za . (^zby^~^ ^<m-l)(*'+l)+;* 

(M' “t~ ^XM' + I + V + l)... (h- + I "1* — 2)(v + l)} 

^ zax^ . (zb . 

< (■’> 

since n+i+(w — 2)(v+i)>i + «?— 2 


m — \. 


But from (10) and (13) we obtain 

if* 1/2 

I (x) — »<,”*> (x) I < I hx'’ w’^ {y) ^ J 

*^1 — s — ) • 

The series S (x) — (x)) and similarly 

7n=0 

the series 2 (^1’"'*'^^ (x) — (x)) are therefore 

«is=0 

absolutely and uniformly convergent in any finite 
interval of x. Hence 

lim (x) = (x), lim (x) = (x). 
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uniformly jy. r. to x. Thus (5) is established. 

• % 

To prove (6) we write 

(x) = 2 {x) - (x)} 2 

+ 2 (x) - (x) } 2. (17) 

n 

Then we have from (13); 

A (x) < 2 2 {x)Y + (x))2] 

+ 2 2 {x) - (x))2 + (x) - (x))]2, 

n 

< 2 (x) f x” ( Ro( j) dy 

•'O v—l 

< 4 r; X" 4b x” ( Ro( j) dj. (18) 

•'0 P—1 

Repeating the inequality (18), we can show by 
induction that 

x> f A I dbx'’+'^ (4^x‘'+i)2 

>’Ro(x)<4ax 

(4b x''+i)’’-i 

(n,+ 1) (n+ 1 + V f i). . . [[i+ 1 +(;>— 2) (v-h i)] 

(19) 


+ ...+ 


} 


< 4^/x'‘ 


(6) is now obvious, because the right side of 
(19) is uniformly bounded in any finite interval 
of X, and is independent of p. 



HYPERBOLIC EQUATIONS II 


Il8 


To prove (7), we must show that for all p 
lim (^) = o (20) 

wi — >• 00 

uniformly w. r. to x. 

Now from (17) we obtain 


•••• 

•>0 

But from (19) we have 

pRo (x) < 4^ A-** 

so that 

pRj (x) ^ lb x” f 4^2 j* dj 

Jo 

< 4a . zb . x” I j« dj 

Jo 

4a . lb . +1 

^ [X -f I 


(21) 


Thus, in general for all i we get 

„Rm M 

4a .{zb AT” +1)” *” 

(H+ i)(fA + I + V + i) ... [g + T -|-"(^•^^Xv+l)] 
4a (zb x” f . 

This proves (20) and consequently (7). 
Proceeding to the limit «? -> 00 in (3) and (4), 
we obtain {ic\ (x), (« — 1,2,..., 00), which 
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are solutions^ cjf (i) and (a), and consequently (x) 
is the solution of (i) § 5.1. 

We shall now prove that (x) is the only 
solution of its kind. If possible, suppose that 

#n C-^) [ with (.x) = is another solution of the 

equation (i) § 5.1., such that 

2 and 2 {v„{x))^ 

n n 

are uniformly convergent. 

Let for all w ^ i ; 

-i- 2 {/'i’"' (.v) - ^„Gv)}2 

n >' 


Then we find 

^ (-v) - "r (-'‘■)} -f 2 (.v) - (x)} 2 

< 2 -I- f 3 J (24) 

It is not difficult to show on the same lines as 
in (22) that 

Urn B„ = o. 


(25) 
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From (7) we know further that Jim = o. 
Consequently, we conclude from (24) that 
2 {»„ (x) - (x)} 2 = o, 2 K (x) — (x)} 2 = o. 

n n 

(26) 

This shows that //„ (x) = (.v) and (x) 

= (x) for every so that the solution (x) 

of (i) §5.1 is unique. 

5-3. The case when g„ (x, x) ^ o. 

In this case the doubly infinite system of non- 
linear integral equation (i) and (3) § 5.1 can be 
written, when (x) is put for z/', (x), 

Cx 

"nW = E.n (^.j) F„{j; K U)} 

K M = gn F\{x-, (a-); Vr W (^) 

+ ^-gn ('•<•> Fn {.r> (j); ( j) ! 

In addition to the hypothesis (2) § 5.1 and 
(8), (9) § 5 . 2, we must now assume that absolute 
constants c,d>o and y, 8 > o exist such that 

2 [gn x)F^{x-,o-,o)\^ < ex'*-, (2) 

2 {al^ (x, x)} ' 

2 {X, X)} ^ • 

n 


(3) 
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We take again the successive approximations 
(x) ^ o, (x) = o, 

(x) = [ (x,j) F„ {j; (j); (j)} dj, 

{x) (x, x) F„ {x; (x); (x)} 

( 4 ) 

We can show that if x lies in a restricted interval 
then for all « > i ; 

lim <”*> (x) = (x); lim (x) -- (x). ... (5) 

m -> 00 m -> 00 

uniformly tv. r. to x, and that 
lim 2; {»!’”’ (x) — //„ (x)}2 = lim 2 {x) 

w->oon m-->-oon 

- (x)} 2 = 0 (6) 

uniformly tv. r. to x. The relation (6) implies that 
(x) and (x) belong to Hilbert space. 

We shall give only a brief outline of the proof: 
the details are as in § 5.2. 

Let 

(x) - 2 (x) - (x)} 2 

n 

+ (x) - V'f (X)} 2]. 


( 7 ) 
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Then from (4) we obtain 

pRm W < r pRm-1 (j) 4y + 2^-^® pRm-1 (^) 

Jo 

+ (* pR„_i (j) dy 

Jo 

< pR„_i (x) + r pR„-i (j) ^ 

Jo 


< (2fl'x« 4- 3/^x’’ I) pR„_i (x), (8) 

where I stands for the integral operator | dj. Thus 

pRm (^) < ( 2 ^.^® + 5 ^-^^’' 1 )” pRo(^) ( 9 ) 


Now we find analogous to (19) that 
,R(, (x) < ^ (x) + (5d'x« + j^x” 1) p_iRo, ... (10) 

< 2 (}dx^ + 5^x’' 1)* ^(x) + (3fllv* tP j^x*" I)*” jRo(x), 
fe =0 

where ^ (.v) = S^zx" 6 cx'^. (ii) 

Similarly 

jRq (x) ^ zfx'^ + 5(7x'‘ (12) 

We assume x to be so small that 

■ (>3) 

It can be shown that on account of (12), the 
last term on the right of (10) tends to zero as -> 00, 
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uniformly w. r. to x provided x satisfies (13), 

and that the sum 2 {'j)dx^ + ')bx'’ I)* <f> (x) remains 

&=0 

bounded. Hence (>;) is uniformly bounded for 
all p. Therefore from (8) we get 

lim ^ o (14) 

W 00 

for all p. 

This shows that for the restricted interval 
(13) the limits 

lim (x) ■-= tt^ (x), lim p^”'^ (x) {x) (15) 

m — > 00 w 00 

exist for every n. We see also that 2 {«„ (x)} ^ 

n 

and 2 (-v)} ^ are convergent. 

n 

These «„ (x) and (x) are solutions of the 
doubly infinite system (i), (2); and it can be proved 
as in §5.2 that the solutions ate unique. 

The system of integral equations 

«» i?<) = r i'n j) (16) 

Jo 

is evidently much easier to deal with than the 
system of integro-differential equations above. It 
can readily be seen that the system (16) has unique 
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solutions belonging to Hilbert space^. 

5*4. Non-linear hyperbolic equation. 

In Chapter II we considered the equation 



P (x, t, u\ 


where P(>r, t, u) was supposed to be an analytic 
function which could be represented in a series 

2 pr {x, f) u''. The method was evidently not 

r=2 

applicable if P contained ^ and ^ also. But the 

system of integro-differential equations studied 
in § 5.3 makes it possible to prove the existence and 
uniqueness of the solution of the equation 



P (x, u, p, q), ... (i) 



d«i 

6/J 


for the boundary conditions 
u (p, t) — u (n, /) = o for all / in o < / ^ T, (2) 


Tor details of § 5 . i — 5 .3, cf. S. M. Sundaram: ‘‘On an 
infinite system of non-linear integro-differential equations.” 
Proc. Ind. Acad. Sc. ^'(^ 939 ) P* 4 ^^^ — 4 ^^* infinite 

system of non-linear integral equations.” Ibid, S (1938) 
p. 238 — 242. 
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and 

u{x, o) = /j (x), ^ (x, o) = /g (x) in o < .V < Jt (5) 

The function e (x) is essentially positive, 
continuous and bounded together with its first and 
second derivatives in o ^ x < Ji : 

o < w < 9 (x) < M, I 

[ Q' (x) I < M, I 9 " (x) I < M. j 

P (x, /, », p, q) is continuous w. r. to all the 
variables and 

P (o, A o, o.,6) = P (.T, /, o, a, o) = o, ( 5 ) 

for all real a and all / > o, so that P can be 
expanded in a series of Sturm-Liouville eigen- 
functions <j>„ (x), [with (/) as in (ii) below] : 

P (x, A u, p,q)^^ P„ [/; iP, (/); jp', (/)} (x), ] 

Pn \ K (0; (0 } ^ a »,p, q)<l>n (^) dx. j 

Jo ) 

We assume further that P is Lipschitzian w. r. 
to u, p, q so that 

{P A P^q) — P A ?)} ^ 

< K, {(« - «)2 + (;> _ ^)2 + (^ _ ^)2} (7) 

where is an absolute constant > o. 

The given functions (x) and /g (x) are also 
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taken to be capable of being expanded in a series 
of <P„ (x); 

/i (x) = S (x), '( 

" i- (8) 

/a (x) = 2 (x), 

n 

We assume further that J[ (x) and /g (x) are 
integrable L^, and write 

2 X „4 + ( 9 ) 

n 

and 

2 (/; o; o) =- j P2 (x, /, o, o, o) dx < Fg (lo) 

• n Jo 

where Fg and Fg are absolute constants. 

For the solution of the equation (i), we write 
u (x, 0 = 2 M'„ (/) <!>„ (x), ... (i i) 

n 

so that the conditions (3) and (8) will be satisfied 
if for all « ^ I we have 

(o) -- ??/ (o) d^. (12) 

Now, if u (x, /) is another function of the form 
(ii), we get 

n 

= [ {P (x, /, », p, q) — P(x, /, u, p, q)} 2 iv 
< K, {£(«-»)Vx+£(p-p)2iv+£(^-^)2ix} 
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< (r‘ + §> 2 

^ ' n n 

< ^^4 {2 - »^„)2 4. 2 « - ^;)2} (13) 

n n 

where K4 — max ( Kj, 

Substituting (6), (ii) in (i), and equating the 
Fourier Coefficients on both sides, we obtain for 
all « > I : 

- -Fn Uy (^); <(01 (14) 

The solution of (14) satisfying the initial con- 
ditions (12) is given by 

cos VX„ / -f -y” sin Vk^ 
v\ V^(/-j)xF„{j;i^,(j);«''(j)}'^, (15) 

(« I, 2, .... co). 

Differentiating (13) jy. r. to / we obtain 
K (0 -- Vx„c„ sin VX„ / + cos V t 

— [‘ cos V X„ (/ -j) X F„ { j; j); K ( j)} (16) 

Jo 

(« = 1, 2, ..., 00). 
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We set 

»« W == K (t) = Kit% (17) 

and get from (15) and (16) the doubly infinite system 

of non-linear integral equations 

«« W = cos V^„ t + sin t 

- sin X F„{j; 

Vn W = - sin yJK t + 4 cos V^„ t 

( t u ( y) 

^ cos X FJj; ;vXy)}dy. 

{n = I, 2, .... co). 

From § 5.3 we see then that one and only one 
solution of the doubly infinite system (18) exists 
for a finite interval of t, such that ul (/) = 2 A.„ 

n n 

(/) and 2 if) — 2 (/)} ^ are uniformly con- 

n n 

vergent, and consequently and are 

integrable L,^ uniformly in t. 

Thus «(x, /) — 2 (/) (•^•) is the required 

n 

solution of (i) satisfying the boundary conditions 
(2) and (3). 

If P {x,t,u,p,q) is defined in the bounded 
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five-dimensional domain 

o < < rt; o < / < T; |»|, \p\, \q\ < N, ... • (19) 

and if P, , P„, P^, P, exist and are bounded, and if 
fi ip^) J" {^) exist and are integrable L.^ so that 
2 4 2 A,„ dl are convergent, then the equa- 

n n 

tion (i) can be solved uniquely in the form (ii) 
so that for sufficiently small values of /, the series 
and 2 {n^nQ)}^ uniformly conver- 

n n 

gent. 

It is not difficult to see that the method can be 
readily extended to the following cases: — 

(a) Instead of the boundary condition (2), 
we can have one of the following boundary con- 
ditions; 




130 


HYPERBOLIC EQUATIONS II 


can have a finite system of equations 

t 



d/2' 


(/ = I, 2, 


= (x, /, ... 

•M^)- (21) 


{c) Instead of two independent variables 
.V, /, we can have more independent variables 
x^, ..., >r„, t so that equation (i) becomes 



= P(.Vi, ...,x^,t,u,p^, ...,p„q) (22) 


where p, = . 

dA-2 

It is obvious that in the extensions (b) and {c) 
the boundary values would have to be generalised 
in a corresponding manner. 
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